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reports . 
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FOREWORD 

Two  studies  completed  during  the  final  year  of  the  grant  are 
reported  herein.  The  first  one  Is  entitled  "Elastostatlc  Fields  In  a 
Crack  Tip  Neighborhood  - Measurements  and  Uniqueness,"  which  is  essential- 
ly the  doctoral  dissertation  submitted  by  J.  H.  Schmidt  to  Texas  A&M 
University.  The  second  study, which  starts  on  page  110,  is  entitled  "On 
the  Mechanical  Behavior  of  Solid  Propellant  under  Transient  Temperatures." 
Both  investigations  deal  in  part  with  the  influence  of  small  flaws  or 
cracks  on  the  mechanical  behavior  of  solid  propellant. 

Work  completed  and  published  during  prior  years  of  the  grant  is 
then  summarized  on  pages  126  - 133.  Abstracts  and  conclusions  are  given 
for  each  of  the  three  studies  covered  on  these  pages. 
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ELASTOSTATIC  FIELDS  IN  A CRACK  TIP  NEIGHBORHOOD  — 
MEASUREMENTS  AND  UNIQUENESS. * 

J.  H.  Schmidt,  Jr. 

ABSTRACT 

Theoretical  and  experimental  investigations  of  stress  and 
displacement  distributions  in  cracked  solids  are  described.  The 
effect  of  damaged  material  in  the  small  highly  stressed  zone  around 
a crack  tip  on  maximum  shear  stress  trajectories  is  first  shown. 
Discussed  next  is  experimental  work  consisting  of  tensile  tests 
conducted  with  long  rectangular  strips  of  solid  propellant  clamped 
along  the  long  edges  and  containing  a central  crack;  displacements 
of  the  strip  were  obtained  through  speckle  diffraction  interfero- 
metry. These  results  are  compared  to  predicted  displacements 
using  a finite  element  technique.  Finally,  a uniqueness  theorem 
for  weak  elastostatic  traction  problems  is  proved;  hypotheses  in 
Kirchoff's  uniqueness  theorem  are  relaxed,  and  uniqueness  is 
proved  for  a variety  of  singular  problems  including  the  crack 
problem. 


♦Based  on  a Ph.D.  Dissertation  submitted  to  Texas  A&M  University, 
May,  1978. 
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INTRODUCTION 

There  are  two  broad  divisions  of  composite  materials,  particulate 
and  fibrous  [1].  This  study  is  concerned  primarily  with  the  appli- 
cability  of  linear  elastic  theory  in  predictinq  elastostatic  states 
near  cracks  in  solid  rccke*  propellant,  a particulate  composite 
material  (Figure  1 ) . 

One's  ability  to  predict  these  states  and  crack  growth  is 
essential  when  predicting  the  service  life  of  particulate  composite 
materials  [2-4], 

The  relevance  of  this  research  becomes  readily  apparent  when 
the  effects  of  large  cracks  in  rocket  motors  are  investigated. 

Cracks  in  rocket  motors  increase  propellant  burning  surfaces,  thereby 
increasing  the  products  of  combustion  with  possible  motor  over 
pressurization  resulting  in  catastrophic  failure.  Additionally, 
propagation  of  the  crack  through  the  propellant  and  motor  exposing 
the  motor  casing  to  hot  gases  can  also  lead  to  failure  [5]. 

A large  number  of  investigators  have  studied  crack  propagation 
in  linear  materials  with  excellent  reviews  of  these  investigations 
given  in  [6,7]. 

Schapery  [2-4]  has  investigated  crack  growth  in  a linear 
viscoelastic  material.  Swanson  [8]  indicated  a reasonable  measure 
of  success  in  the  use  of  Schapery's  theory  in  predicting  both 
initiation  and  subsequent  velocity  of  crack  propagation  in  PRAN 
solid  propellant. 

The  format  of  this  dissertation  follows  the  style  of  the 
International  Journal  o f Frac ture . 
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Previous  investigators  [9-11]  have  documented  the  nonlinear 
behavior  of  solid  propellant;  however,  little  or  no  work  has  been 
done  to  explain  why,  in  some  cases,  a globally  linear  theory  of 
fracture  apparently  predicts  initiation  and  velocity  of  crack  propa- 
gation in  a globally  nonlinear  material. 

The  theories  of  fracture  mechanics  can  be  utilized  to  investi- 
gate, in  part,  this  phenomenon.  The  work  presented  here  begins  with 
the  classical  works  of  Barenblatt  [12]  and  Williams  [13].  Displace- 
ments are  then  obtained,  through  a finite  element  code,  by  solving 
the  classical  linear  elastic  equations  applied  to  a cracked  biaxial 
strip  — a long  rectangular  strip  clamped  along  the  long  edges  and 
containing  a central  crack.  These  numerically  obtained  displacements 
are  compared  to  displacements  measured  on  a cracked  biaxial  strip 
of  PBAN  propellant  through  a relatively  new  laser  scheme,  speckle 
diffraction  interferometry. 

Questions  about  uniqueness  of  solutions  to  elastostatic  traction 
boundary  value  problems  are  addressed.  A uniqueness  theorem  for 
bounded,  locally  star-shaped,  regular  domains  is  proven.  The 
boundary  value  problem  is  cast  in  a "weak"  sense  (to  be  explained 
in  a following  chapter)  requiring  continuous  displacements  and 
absolutely  integrable  (in  the  Lebesgue  sense)  tractions  over  the 
boundary. 
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A REVIEW  OF  LINEAR  ELASTIC  FRACTURE  MECHANICS 
General  Development 

Since  the  beginning  of  the  twentieth  century  there  has  been  an 
increase  in  the  number  of  articles  published  on  fracture.  There  are 
three  basic  types  of  stress  fields  near  crack  tips  [14],  each  asso- 
ciated with  a specific  mode  of  deformation,  as  illustrated  in  Figure 
2.  It  is  Mode  I that  we  are  concerned  with  in  this  work.  For 

completeness,  however,  the  stresses  (o  , o , t ) and  displacements 

x y *y 

(u,  v,  w)  for  all  three  modes  are  listed  below  as  given  by  Irwin  [15] 
(with  r,e  the  local  polar  coordinates  in  the  x-y  plane). 


MODE  I 
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The  parameters  Kj,  Kjj,  and  Kjjj  are  the  stress  intensity  factors 
corresponding  to  the  three  modes  of  deformation  while  v is  Poisson's 
ratio  and  v is  the  modulus  of  rigidity  (shear  modulus).  It  is 
important  to  note  that  while  the  stress  intensity  factors  are  inde- 
pendent of  the  coordinates  r and  0,  they  are  proportional  to  the 
loading  forces  and  also  depend  on  the  configuration  of  the  body  [14]. 

There  are  two  physically  distinct  types  of  plane  elastostatic 
problems  that  we  shall  be  concerned  with;  plane  stress  and  plane 
strain  (plane  deformation). 

The  underlying  assumption  of  the  plane  stress  problem  is  that  the 
surfaces  of  the  object  are  stress  free.  This  assumption  together 
with  the  equilibrium  equation  implies  — rr  = 0 along  the  strip  surface 
(where  the  z-axis  is  normal  to  the  surface).  This  suggests  that,  for 
a thin  strip  (t/B  <<  1 where  t = strip  thickness  and  B = any  in-plane 
geometric  dimension),  oz  = 0 [16]. 

The  plane  strain  or  plane  deformation  problem  is  characterized 
by  u = u(x,y),  v = v(x,y)  and  w = 0 where  u,  v,  and  w represent 
the  displacements  [16]. 

Displacements  given  by  equations  (1),  (2),  and  (3)  are  for  the 
case  of  plain  strain  and  can  be  regarded  as  a good  approximation  in 
the  region  defined  through  r/B  <<  1 where  6 is  the  distance  from  the 
crack  tip  to  the  nearest  geometric  feature  [17]. 

In  the  cracked  plane  problem,  the  stress  field  away  from  the 
crack  tip  can  often  be  considered  to  be  in  a state  of  plane  stress, 

| while  the  stress  field  in  the  crack  tip  neighborhood  suggests  a state 

of  plane  strain  [4,  17].  The  plane  stress  and  plane  strain  solutions 
differ  only  in  the  displacement  fields;  i.e.,  the  in-plane  stress 
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fields  are  identical  [18].  In  obtaining  equations  (1) — (3),  Irwin 
[15,  20-22]  used  complex  stress  functions  Z and  V - based  on  the 
method  of  Westergaard  [19]  that  satisfy 

ox  = ReZ  - ylmZ' 

oy  = ReZ  + ylmZ'  (4) 

'xy  = - yReZ'  • 

Westergaard  [19],  using  complex  stress  functions  (for  a modern 
treatise  of  complex  stress  functions,  see  for  example,  Timoshenko 
and  Goodier  [18],  who  closely  follow  Muskhelishvili  [23]),  solved 
the  two-dimensional  problem  of  a finite  through-crack  in  an  infinite 
sheet  with  ox  = o at  |z|  = «°. 

It  was  shown  by  Sih  [24]  and  by  Eftis  and  Liebowitz  [25]  that 

the  formulation  by  Westergaard,  which  was  subsequently  used  by  Irwin, 

was  invalid  for  a class  of  problems  dealing  with  an  infinite  medium 

containing  finite  cracks  with  loads  applied  at  infinity.  Both  Sih, 

and  Eftis  and  Liebowitz,  proved  that  both  o and  a (equation  (4)) 

x y 

must  be  modified  by  an  additive  constant.  For  example,  in  the 
opening  mode  problem,  with  the  following  boundary  conditions  applied 
to  an  infinite  sheet  containing  a stress  free  interior  crack 
[25], 
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the  following  approximations  for  the  plane  strain  crack  tip  stress  and 
displacement  fields  are  obtained  [25] 

°X  * (7n7)T7?  C°S  ? " (S<M  ^)(Sln  TJ]  ‘ 0 ' k)<’ 

°y*^T7?  cosl  ['  * (s1"|)(s,nf)] 


1/-j..  0\/ 0\/ 3o, 


°xy  " 7<sin  p(cos  p(cos 


HU  a K, 


[*1 


cos  -0  (l  - 2v>)  + si 


. 2 o I 
i n w 


- ?•  (1  - v^)(l  - k ) r cos  o 


“ Ki  [h]  ' sin  I [(2  - 2v)  - cos2 1] 


+ ^v(l  +v)(l  - k)  r sin  o , 


For  the  equal  biaxial  stress  field,  k = 1,  it  is  clear  that  the 
additive  constant  tern  in  the  ox,  u,  and  v equations  (equation  (6)) 
are  zero. 

Inglis  [26],  also  in  the  context  of  classical  elasticity  theory 
presented  the  solution  of  an  infinite  sheet  containing  an  elliptical 
cavity  with  the  sheet  subjected  to  a uniform  stress  field  applied  at 
infinity  normal  to  the  ellipse  major  axis.  Inglis  [26]  then  obtained 
the  solution  to  the  crack  problem  by  requiring  the  minor  axis  of 
the  ellipse  to  approach  zero.  Muskhel ishvi 1 i [23]  also  in  the 
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framework  of  classical  elasticity— sol ved,  among  others,  the  two- 
dimensional  problem  of  a semi-infinite  straight  cut  in  an  infinite 
sheet  with  the  faces  of  the  cut  subjected  to  internal  pressure. 

A.  A.  Griffith  [27],  in  a fundamental  paper  on  the  theory  of 
cracks  in  brittle  fracture,  treated  the  problem  of  a finite  crack  in 
an  infinite  sheet  acted  upon  by  a uniform  tensile  stress,  o,  at 
infinity.  A critical  stress,  oQ,  was  determined  at  which  the  crack 
began  to  expand.  Griffith's  results  were  based  on  the  work  of 
Inglis  [26]  and  an  extension  of  the  "theorem  of  minimum  potential 
energy"  by  taking  into  account  the  increase  in  surface  energy  which 
occurs  during  crack  growth. 

Williams  [13]  studied  the  elastic  stress  distribution  in  the 
neighborhood  of  a crack  under  arbitrary  loadings  and  showed  the 
stresses  to  be  of  order  0(r"^^)  as  r -+■  0.  Williams'  solution 
was  based  on  his  earlier  work  [28]  dealing  with  singularities  in 
angular  corners.  In  [28],  Williams,  using  an  Airy's  stress  function 
and  a displacement  function  in  the  sense  of  Coker  and  Filon  [29], 
solved  the  angular  corner  problem  with  stress  free  surfaces.  In  his 
later  work,  Williams  [13]  allowed  the  included  angle  of  the  angular 
corner  to  approach  2n.  The  resulting  stress  and  displacement  fields 
for  the  opening  mode  (Mode  I)  are  given  for  plane  strain  in  polar 
coordinates  by  the  following  expressions.  (Williams  solved  the 
plane  stress  problem;  however,  as  noted  earlier,  in  the  crack  tip 
neighborhood  a condition  of  plane  strain  exists.) 
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where  the  quantity  6 represents  the  distance  from  the  crack  tip  to  the 
nearest  geometric  feature.  It  should  be  emphasized  that  Williams1 
solution  (equation  (7))  yields  a square  root  singularity  in  stress 
regardless  of  external  loadings  and  shape  of  the  crack  and  body. 

Williams'  solution  (equation  (7))  and  the  corrected  Irwin  solu- 
tion (equation  (6))  are  equivalent  in  the  crack  tip  neighborhood  for 

symmetric  loadings  when  the  constants  a.  and  K.  are  related  through 

-K  * 

a,  = — - and  a0  = - Similarly,  Williams'  antisymmetric 

1 2 A 

stress  case  (not  shown  here)  is  equivalent  to  Irwin's  Mode  II  in  the 

Kj  i 

crack  tip  neighborhood  under  the  substitution  b.  = . These 

✓ST 

equivalences  are  shown  in  Appendix  A for  reference  purposes. 

Although  Williams'  expansion  is  not  suitable  for  internal  cracks, 
the  method  employed  by  Williams  is  readily  extended  to  this  casr  when 
the  complex  variable  form  is  used  [30,  31]. 


1? 

Limitations  of  the  Linear  Elastic  Solution 

In  these  early  investigations  the  equations  of  linear  elasticity 
were  solved;  i.e.,  Hookes'  law  and  small  displacement  qradient  theory 
were  assumed.  In  the  linear  theory,  particles  close  together  relative 
to  the  size  of  geometric  features  are  assumed  to  remain  close  together 
after  deformation;  however,  when  cracking  occurs,  this  close  neigh- 
borhood relation  fails  to  hold.  The  solution  to  the  linear  elastic 
problem  in  these  early  investigations  led  to  strains  and  stresses 
which  were  no  longer  small  compared  to  unity;  indeed,  the  strains  and 
stresses  become  unbounded  as  the  crack  tip  is  approached,  and  the 
solution  is  simply  not  valid  in  this  region  [32]. 

There  are  three  conditions  which  lead  to  singular  solutions  in 
an  elastic  medium  [33];  geometric  discontinuities,  material  dis- 
continuities, and  load  discontinuities.  Geometric  discontinuities 
are  of  principal  concern  in  this  work  and  for  discussions  of  this 
type  of  singularity  the  reader  is  referred  to  the  work  of  Sternberg 
[33],  Sternberg  and  Koiter  [34],  and  England  and  Spencer  [35]. 

The  application  of  linear  elastic  solutions  to  the  crack  problem 
is  validated  by  the  exclusion  of  a small  zone  of  such  size  and  shape 
that  the  remaining  region  satisfies  linearity.  This  zone  is  commonly 
referred  to  in  the  literature  as  the  plastic,  failure,  process,  or 
damage  zone.  This  dissertation  will  refer  to  this  zone,  when  in  the 
form  of  a thin  layer,  as  the  failure  zone  of  length  o as  shown  in 
Figure  3. 

A criterion  (e.g.  Griffith's  energy  criterion  [26])  is  needed  to 


determine  when  crack  growth  occurs;  i.e.  when  the  crack  boundary 
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undergoes  a change.  In  this  work  we  will  be  interested  in  stationary 
cracks. 

To  date,  the  modelling  of  crack  tip  material  has  been  restricted 
to  two  approaches:  a)  the  "particle  approach"  where  individual  atoms 
and  interatomic  forces  are  idealized  through  some  mathematical  model, 
and  b)  the  "continuum  approach"  where  mass  is  not  concentrated  in 
particles,  but  rather  in  small  volumes.  This  latter  model  filters  out 
many  of  the  complexities  of  the  former  approach;  however,  in  dealing 
with  cracks,  the  microstructure  behavior  is  helpful  in  the  crack  tip 
neighborhood  in  the  interpretation  of  stress  behavior  in  the  failure 
zone.  It  is  in  this  crack  tip  neighborhood  that  the  continuum  model 
falls  short.  Barenblatt  [12]  combines  these  two  models. 

We  shall  begin  by  following  the  lead  of  Barenblatt  [12],  Dugdale 
[36],  and  others,  and  postulate  that  the  stresses  in  the  elastic 
medium  must  remain  finite— a physical  condition.  Following  Schapery 
[2]  we  designate  the  stress  in  the  failure  zone  as  the  failure  stress 
distribution,  o^.. 

Schapery  [2-4]  and  Goodier  [32]  proposed  a crack  tip  model 
mathematically  similar  but  conceptually  different— in  the  interpreta- 
tion of  the  finite  stress  condition— to  Barenblatt's  crack  tip  model. 
Dugdale  [36],  in  a study  of  the  plastic  zone  at  the  end  of  a finite 
slit  in  an  ideal  elastic-plastic  material,  used  the  finite  stress 
postulate  in  obtaining  the  plastic  (failure)  zone  length. 

It  is  interesting  to  note  that  Schapery  [2]  proved  the  equiva- 
lence of  the  cohesive  force  model  and  the  elastic  model  of  Williams 
[13]  when  the  "failure  zone"  length  of  the  cohesive  force  model 
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approaches  zero.  Rice  [37]  indicated  that  the  energy  criterion  for 
fracture  of  the  cohesive  model  reduced  to  the  Griffith  theory, 
regardless  of  the  force  attraction  law,  whenever  the  "failure  zone" 
was  small  compared  to  characteristic  dimensions. 

Referring  to  Figure  3,  the  "apparent  crack  tip",  following 
Schapery  [2],  is  located  by  definition,  at  the  point  x = -a,  (though 
this  location  is  often  called  the  crack  tip  in  the  literature),  while 
the  crack  tip  is  located  at  x = 0.  All  theoretical  crack  equations 
referred  to  in  this  work  are  written  with  respect  to  a coordinate 
system  located  at  the  crack  tip;  however,  because  of  the  difficulty 
in  experimentally  locating  the  crack  tip,  all  experimental  data  is 
recorded  relative  to  the  apparent  crack  tip  (cf.  the  chapter  comparing 
theory  and  experiment).  In  order  to  study  the  effect  of  the  failure 
zone  on  the  stress  and  displacement  fields  we  shall  follow  Schapery 
[2],  who  removed  the  stress  singularity  by  considering  a failure  zone 
ahead  of  the  apparent  crack  tip.  The  model  used  was  obtained  by 
superposing  the  solutions  of  Williams  [13]  and  Muskhel ishvi 1 i [23]. 

Schapery  further  extended  the  elastic  solution  to  that  of  a 
viscoelastic  medium  via  the  correspondence  principle.  Although  the 
undamaged  portion  of  the  continuum  was  assumed  to  be  a linear, 
isotropic,  viscoelastic  medium,  the  material  within  the  zone  could  be 
highly  nonlinear;  indeed,  it  need  not  even  be  a continuum,  in  that  it 
could  consist  of  strands  of  material  as  in  Figure  3. 

In  Figure  4 the  reaction  of  the  failure  zone  on  the  linear 
continuum  is  shown,  and  is  designated  the  "failure  stress",  a^. 
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AN  ANALYSIS  OF  THE  FINITE  STRESS  CRACK  PROBLEM 
General  Discussion 


A detailed  analysis  of  the  combined  Barenblatt  and  Williams 
problem  will  now  be  presented  so  that  the  effects  of  the  failure  zone 
on  the  stress  and  displacement  fields  can  be  determined. 

Since  closed  form  solutions  to  the  crack  problem  only  exist  for 
the  case  of  an  infinite  body,  we  will  restrict  our  study  to  this  case. 
In  order  to  analyze  the  finite  stress  problem  we  shall  first  consider 
the  stress  and  displacement  fields  of  two  singular  problems.  The 
first  problem,  solved  by  Barenblatt,  involves  an  infinite  body  with  a 
semi-infinite  cut  that  is  symmetrically  loaded  with  normal  stress, 
along  the  surface  of  the  cut.  The  second  one  is  the  elastic 
solution  to  the  crack  problem,  given  by  Williams  [13]. 

The  Barenblatt  Problem 


From  Barenblatt  [12]  we  obtain  the  solution  to  the  first.  The 
solution  is  essentially  that  given  by  Muskhelishvili  (231  restricted 
to  normal  stresses  acting  along  the  surface  of  the  cut. 


= 4 Re  $(z) 


°y  " ioxy  = + + <z  " 


(8) 


with  displacements 


2w(u'  + i v ' ) = k $(z)  - u>(z)  - (z  - z)  <t> ( z ) 


(9) 
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where  for  plane  strain,  k = 4 - 3v. 

The  analytic  functions  w,  $,  and  w are  expressed  through 


«(z)  = i.'(z)  = = w'(z)  = 


1 

2.i  i/z 


/Vt  of(t) 

y v- dt 


♦ (z)  = co(z)  = -X  f 0f(t)  In 

0,1  J T /t  - /z 


dt 
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The  following  failure  stress  distributions,  cy,  will  be 
considered. 


a)  of(t)  = tc [0 ,«) 


(Ha) 


= 0 otherwise 


b)  of (t) 


T 1 * "m 


(lib) 


= 0 otherwise 


These  stress  distributions  are  for  a crack  along  the  positive 
real  axis  since  Barenblatt's  [11]  and  Muskhel ishvil i 's  [23]  work 
are  cast  in  this  way.  The  solution  of  Williams  [13]  is  for  the 
crack  along  the  negative  real  axis 

For  a crack  along  the  positive  real  axis,  the  branch  cut  for  the 

1/2 

"multivalued"  functions  In  (•)  and  (•)  ' is  also  along  the  positive 
real  axis.  This  poses  no  real  problem,  even  though  we  need  to  super- 
pose the  two  solutions.  We  can  define  z 3 -x  + iy  which  in  effect 
places  the  crack,  and  thereby  the  branch  cut,  along  the  negative  real 


axis.  Additionally,  by  physical  reasoning,  the  shear  stresses  and 
the  displacements  in  the  x-direction  of  the  Barenblatt  solution  must 
undergo  a sign  change. 

We  seek  the  elastostatic  state  corresponding  to  the  failure 
stress  distribution  of  equation  (11a).  We  need  to  find  the  analytic 
function,  <j> ( 2 ) . Towards  this  end  consider 


(16) 
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Equations  (8),  (9),  (14),  (15),  and  (16)  qive  the  state  of  stress 
for  the  Barenblatt  problem  when  is  defined  as  in  equation  (11a). 

Now  we  turn  our  attention  to  the  Barenblatt  problem  when  the 
failure  stress  distribution  is  given  by  equation  (lib). 

As  before, we  must  find  <jv  for  the  given  failure  stress  distribu- 
tion. Consider  the  following. 


jr  ' (t  - z)ln(i5-L£) 

Si  - Si 


2 _2, 


+ 2/zVt 

St  - Si 


- Sill.  . ills,  L (,  . 

a a ( a Si  - Si 


whenever  z + t.  Clearly 


Al  - ^)ln(^--ii)dt  = « [(1  - J)in(1---^:) 
Jo  a Si  - Si  L 1 - SiTa 


- (I^Za))  in fJ-LjSi)  + | Siu  - t)V2 
2 1 - Sz/n  3 " 


whenever  z f a.  The  resulting  function  is 


■ £ ["  - £>*  * f 
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which  is  analytic  everywhere  in  the  complex  plane  except  along  the 


branch  cut  (the  positive  real  axis).  Because  of  the  analyticity  of 
$ • 
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>"(*)  = [ ln(LjL2^)  - 2(-==)  - |( Vf} 

2nl  « L 1 - /zjTa  fi/a  3 (z/a)3/2 
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Equations  (8),  (9),  (19),  (20)  and  (21)  completely  define  the 
stress  and  displacement  fields  to  the  Barenblatt  problem  with  a failure 
stress  distribution  given  as  the  ramp  function  in  equation  (lib). 


The  Williams  Problem 


The  polar  stresses  and  displacements  (equation  (7))  given  by 
Williams  must  now  be  written  as  rectangular  components  of  stress 
and  displacement  so  that  they  may  be  superposed  with  the  rectangular 
stress  and  displacement  components  given  in  the  Barenblatt 
solution. 

We  seek  u and  v,  the  rectangular  components  of  the  Williams 
solution  (equation  (7)).  By  a coordinate  transformation  we  obtain 

u = (cos  o)ur  - (sin  o)u0 

(22) 

v = (sin  o)u  - (cos  o)u . 

r u 

Timoshenko  and  Goodier  1181  related  the  rectangular  stresses, 
(v  °y‘  \y).  t0  the  polar  Stresses  (or,  o0»  xr0)  through 
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2 2 

ox  = or(cos  o)  + o()(sin  o)  - 2y0  sin  o cos  o 

p p 

0y  = or(sin  o)  + oQ( cos  o)  + 2Tr0  sin  o cos  0 (23) 

2 2 

t = (o  - o„)sin  0 cos  0 + t „(cos  0 - sin  0) 
xy  r 0 r0'  ' 


by  equilibrium  considerations.  Equation  (23)  can  also  be  obtained 
through  a coordinate  transformation  of  the  second  order  stress 


tensor. 


The  Combined  Problem 

Next,  we  add  the  solution  of  Williams  [13]  (equations  (7),  (22), 
(23)),  for  the  case  of  equal  biaxial  loading  , i.e.  = 0 in  equation 
(7),  and  Barenblatt  [12]  (equations  (8),  (9),  (10)).  The  stress  state 
obtained  from  the  addition  of  two  stress  states  is  an  application  of 
the  method  of  superposition.  Superposition  is  possible  as  a result 
of  the  linearity  of  the  system  of  governing  partial  differential 
equations,  cf.  Sokolnikoff  [16].  This  addition,  together  with  the 
finite  stress  postulate,  yields  a relatively  simple  expression  re- 
lating Kj  and  of  [4,32], 


For  the  failure  stresses  considered  in  this  work  (equations  (11)), 
the  relationship  between  Kj  and  cy  (equation  (24))  results  in 

K.  • 2(2V/2o  <»*> 

I ' n ' m 

for  the  failure  stress  in  equation  (11a),  and 
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KI  = 3 <¥*  °m  (25b> 

for  the  failure  stress  associated  with  equation  (lib). 

The  elastostatic  state  corresponding  to  the  stress  distributions 
of  equations  (11)  is  obtained  through  the  superposition  of  equations 
(7),  (8),  (9)  (where  a^  = - Kj/2h  in  equation  (7))  together  with  the 
finite  stress  conditions  of  equations  (25a)  and  (25b).  Figures  5 and 
6 compare  the  grid  deformations  for  these  combined  problems,  corres- 
ponding to  equations  (11a)  and  (lib)  respectively,  to  the  grid 
deformations  given  by  the  Williams'  solution.  Particular  attention 
is  drawn  to  the  deformations  near  the  crack  tip,  which  show  the 
opposite  faces  of  the  crack  surface  closing  smoothly  in  a cusp-like 
fashion,  a result  due  to  the  finite  stress  condition  (13).  The 
nondimensional ized  normal  stresses  acting  along  the  crack  plane  are 

shown  in  Figure  7,  which  graphically  illustrates  the  continuity  of 
°v 

— at  the  crack  tip. 

°m 

Since  the  stress  intensity  factor  is  determined  by  geometry 
(including  crack  size)  and  loading  conditions  [14],  it  is  reasonable 
to  normalize  the  stress  fields  corresponding  to  the  different  failure 
stress  distributions  (equations  (11a)  and  (lib))  with  respect  to  the 
same  stress  intensity  factor.  This  is  accomplished  for  the  ramp  and 
constant  failure  stress  distributions  by  equating  equations  (24)  and 
(25)  and  by  associating  o with  a specified  maximum  strength. 

This  results  in  different  failure  zone  lengths  which  are 

related  through  a I = v a I . The  effect  of 

' ramp  I constant 

the  two  failure  stress  distributions  are  displayed 


id  grid  deformations  of  the  combin 
ams'  problem  for  a ramp  failure  s 
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in  Figure  8, which  indicates  the  dominance  of  the  Williams  singularity 
term  along  the  plane  of  the  crack  after  a distance  of  approximately 
5a. 

Figures  9-12  compare  the  isochromatic  fringe  patterns*  of 
the  singular  Williams  solution  to  that  of  the  combined  problem  for 
both  the  constant  and  ramp  failure  zone  stresses.  The  tilting  of  the 
isochromatic  lines  is  due  solely  to  the  cohesive  zone  and  cannot 
theoretically  be  answered  through  the  addition  of  a constant  term  to 
the  singular  term  in  Williams'  c>x  expansion,  since,  as  previously 
discussed,  this  term  is  zero  for  equal  biaxial  loadings. 

For  the  case  of  uniaxial  tension  applied  to  a plate  with  an  edge 
crack, the  tilting  effect  of  the  isochromatic  fringes  was  explained  by 
previous  investigators  [38-40]  through  the  addition  of  a constant  to 
the  singular  term  of  the  Williams'  expansion.  That  the  constant  term 
in  the  expansion  must  be  included  is  a theoretical  reality  [24, 

25];  however,  whether  or  not  the  isochromatic  fringe  tilting  is  due 
exclusively  to  this  term  with  no  regard  to  the  cohesive  zone  depends 
on  the  scale  of  this  zone. 

In  Figures  9 and  10  the  isochromatic  fringes  touch  the  crack 

plane  at  x/a  = -1.0  because  of  the  finite  stress  discontinuity 
o o 

(-■*--*.0  as  x/a  -*•  -1.0  from  the  left  and — *■  ■+  1 as  x/a  * -1.0 

O O _ 

n m 

from  the  right)  of  the  constant  failure  stress  distribution.  These 

illustrations  show  that  the  shape  of  the  isochromatic  fringes  in  the 

crack  tip  neighborhood  is  very  dependent  on  the  failure  stress 

distribution. 

* 

These  patterns  are  the  maximum  shear  stress  trajectories. 


Williams  Singular  Solution 


1 
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FINITE  ELEMENT  SOLUTIONS  OF  THE  BIAXIAL  STRIP 
Motivation 

The  most  crucial  part  of  this  work  is  the  experimental  verifi- 
cation of  analytical  results.  Inasmuch  as  linear  theory  is  used 
we  come  to  the  question  of  its  validity  in  the  practicing  engineering 
environment. 

The  cracked  infinite  sheet  problem  could,  in  principle,  be 
approximated  in  the  laboratory.  This  approach  would  be,  at  best, 
cumbersome.  An  alternative  to  this  approach  is  the  numerical  solu- 
tion to  a more  easily  handled  geometry  with  more  readily  applied 
boundary  conditions. 

The  exact  solution  to  many  boundary  value  problems  of  elasticity 
in  engineering  often  presents  formidable  mathematical  difficulties, 
due  to  the  complicated  form  of  the  physical  boundary  conditions. 

It  is  noteworthy,  in  technical  applications  of  the  theory  of  elasti- 
city, that  for  many  problems  only  a mathematical  approximation  to  the 
physical  boundary  conditions  exist.  The  resulting  solution  then 
represents  only  an  approximation  to  the  physical  problem.  The  finite 
element  method  can  provide  solutions  based  on  an  improved  approxima- 
tion of  the  physical  boundary  conditions. 

Finite  Element  Analysis 

Numerically  obtained  displacement  fields  and  grid  deformations 


were  compared  to  experimentally  determined  displacement  fields 
and  grid  deformations.  These  numerical  solutions  were  calculated 
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by  means  of  an  existing  finite  element  code. 

Due  to  the  symmetry  of  the  chosen  test  geometry  only  a quarter 
of  the  biaxial  strip  was  analyzed.  Following  Conrad  [41],  The  Texas 
Grain  Analysis  Program  (TEXGAP)  [4?]  (a  finite  element  program 
originally  developed  for  analyzing  solid  rocket  propellant),  was 
used  to  numerically  analyze  the  quarter  biaxial  strip.  There  are 
standard  texts  which  discuss  the  finite  element  method  (cf.  Prenter 
[43],  Strang  and  Fix  [44],  Oden  and  Reddy  [45]).  For  many  boundary 
value  problems,  an  energy  inner  product  exists  (related  to  the 
potential  energy  of  an  appropriate  physical  system)  whose  minimum 
with  respect  to  some  admissible  class  of  functions  is  the  solution 
to  the  boundary  value  problem.  This  energy  formulation,  when  applied 
to  a body  which  is  approximately  using  discrete  volume  or  area 

elements,  is  referred  to  as  the  finite  element  method. 

Three  numerical  problems  were  solved  for  a cracked  biaxial  strip 
with  dimensions  of  2 inches  by  7 inches.  A comparison  of  each  of  the 
three  numerical  problems  with  experimental  data  is  recorded  and 
discussed  following  a discussion  of  the  experimental  work. 

First,  the  problem  of  a 1.5  inch  crack  in  a biaxial  strip  with 
no  failure  zone  (the  square  root  singular  solution)  was  solved.  The 
second  type  of  problem  investigated  was  a superposition  of  the  solu- 
tion to  the  first  problem  (with  crack  lengths  of  1.7  inches  and  1.9 
inches),  and  the  solution  to  the  Barenblatt  infinite  sheet  problems 
discussed  previously  (with  failure  zone  lengths  of  .1  inches  and  .2 
inches  respectively).  At  this  point  it  was  felt  that  a more  refined 
numerical  problem  not  requiring  the  superposition  of  a geometrically 
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finite  problem  with  that  of  a geometrical ly  infinite  problem  should 
be  solved.  Therefore,  the  third  problem  was  solved  by  linking  the 
failure  stress  distribution  with  the  TFXGAP  solution  through  the 
finite  stress  condition  (equation  (24))  as  shown  in  Appendix  B. 

The  plane  stress  mode  of  the  TEXGAP  code  was  used  since,  as 
previously  indicated,  the  sheet,  except  for  a small  neighborhood 
about  the  crack  tip,  is  approximated  by  the  plane  stress  hypothesis. 

The  generation  of  the  crack  elements  by  the  TEXGAP  program  is 
discussed  by  Conrad  [41]  and  also  by  Dunham  and  Becker  [42].  The 
elements  generated  by  TEXGAP  in  solving  the  quarter  biaxial  strip 
problem  in  this  work  are  rectangular  and  are  shown  later  when  a 
comparison  of  theory  and  experiment  is  made. 


36 


FRACTURE  TESTS 
Samples  and  Equipment 

Composite  solid  propellant  consists  of  a rubber  matrix  filled 
with  hard  particles.  The  propellant  used  in  this  investigation  was 
polybutadiene  acrylonitrile  (PBAN)  containing  66%  by  volume  ammonium 
perchlorate  with  a particle  size  range  of  20-100  microns  and  10%  by 
volume  aluminum  powder  with  particle  size  on  the  order  of  10  microns. 
The  samples  were  provided  by  the  Air  Force  Rocket  Propulsion  Laboratory. 
They  were  stored  in  a controlled  climatic  environment  of  75°F  and 
25%  relative  humidity  prior  to  testing. 

The  geometry  of  the  test  specimen  chosen  for  the  experiments 
reported  herein  was  that  of  the  biaxial  strip  as  shown  in  Figure  13. 

The  nominal  sizes  of  the  specimen  were  2 inches  X 7 inches  with 
thicknesses  of  .125  and  .2  inches. 

Samples  were  bonded  between  two  stainless  steel  bars.  The 
bonding  jig  used  to  position  the  sample  relative  to  the  steel 
bars  is  shown  in  Figure  14.  The  jig  was  designed  so  that  a micro- 
meter could  be  used  to  center  the  sample  to  within  ^.002  of  an  inch. 

One  face  of  the  sample  was  bonded  to  one  bar  and  allowed  to 
cure;  then,  the  sample  was  removed  from  the  jig,  turned  up  side  down, 
and  the  opposite  face  bonded  to  the  remaining  stainless  steel  bar. 

The  two  steel  bars  were  kept  from  compressing  the  sample  during  this 
second  bonding  operation  by  two  precision  1/2  inch  X 2 inch  pins 
placed  between  the  two  opposing  steel  bars.  The  epoxy  used  through- 
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out  this  work  was  Versamid  125  hardener  and  R-815  resin  with  a hard- 
ener to  resin  mixing  ratio  of  2 to  3 (by  volume).  All  specimens  had 
a curing  time  of  over  24  hours.  Prior  to  bonding, the  stainless  bars 
were  cleaned  with  acetone  and  then  treated  with  M-Prep  Conditioner 
A and  B. 

After  samples  were  bonded, a rigid  aluminum  plate  was  attached 
to  the  steel  bars  preventing  premature  sample  straining. 

Flat  white  Krylon  paint  was  sprayed  on  the  specimen  surface 
prior  to  testing  to  enhance  interferometric  fringes— the  experimental 
method  employed  to  measure  displacements.  During  the  painting 
process  the  paint  can  was  held  approximately  30  inches  from  the 
specimen  surface  in  an  attempt  to  minimize  +he  paint  thickness  by 
mist  spraying  the  specimen.  The  effect  of  the  painting  technique  is 
addressed  in  a later  section. 

Finally,  a crack  was  cut  using  a vertical  milling  machine  to 
press  a razor  blade  through  the  sample.  A crack  cut  in  this  fashion 
could  be  placed  within  .005  of  an  inch  from  the  specimen  center. 

Displacement  Measurement  Methods 

We  now  address  the  experimental  problem  of  displacement  field 
measurements.  Of  the  methods  used  for  measuring  displacements,  the 
following  were  considered  most  promising:  (i)  Moire',  (ii)  a 
pattern  recognition  technique  utilizing  a microdensitometer,  and 
(iii)  speckle  diffraction  interferometry. 

A close  examination  of  the  Moire'  technique  reveals  some  obvious 
disadvantages,  the  foremost  being  the  need  to  apply  a grid  to  the 
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specimen  surface.  The  grid  pattern  must  be  of  the  highest  quality 
with  regard  to  line  straightness,  uniformity  of  pitch,  contrast,  and 
uniformity  of  line  to  space  ratio.  At  the  same  time, the  pattern 
must  have  the  necessary  line  density  to  achieve  required  sensitivity 
146]. 

There  are  several  methods  available  for  grid  application,  among 
them  are  photographic  etching,  machine  ruling,  and  transferable  grids. 

Photoetching  is  a delicate  procedure, and  imperfect  photographic 
reproduction  of  the  master  grid  can  produce  large  errors  in  strain. 
Moreover,  sensitizers,  developing  materials,  and  dyes  used  in  the 
photoetching  procedure  (e.g.  ammonia  perchlorate)  act  as  desiccants, 
altering  the  chemical  and/or  physical  properties  of  the  sample. 
However,  Hart  [47]  employed  a method  which  for  propellant  samples 
apparently  overcomes  these  difficulties  by  using  a diazo-type  photo- 
sensitive material  together  with  ammonia  vapors.  Hart,  using  a 
grid  of  500  lines/in,  found  that  an  average  strain  of  9.2% 
was  required  to  obtain  meaningful  fringes.  By  comparison,  using 
speckle  diffraction  interferometry, the  present  tests  had  average 
strains  of  .2%. 

The  most  promising  grid  application  is  that  of  grid  transfer. 

The  best  available  grid  transfers  have  grid  densities  on  the  order 
of  2000  1 ines/in.  [48]  , resulting  in  a displacement  sensitivity  of 
500  pinches.  On  the  other  hand,  the  interferometric  technique  has  a 
displacement  sensitivity  of  40  uinches.  Additionally,  excessive 
specimen  distortion  can  occur  due  to  adhesion  between  the  transfer 
grid  and  specimen  when  the  adhesion  modulus  is  on  the  order 
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or  greater  than  the  specimen  modulus. 

The  microdensitometer  technique  at  first  seems  very  promising; 
however.it  is  very  expensive  both  in  initial  expense  and  subsequent 
computer  expense.  Furthermore,  it  does  not  appear  to  be  as  accurate 
as  the  interferometry  technique.  Essentially,  two  film  negatives 
are  needed,  one  before  and  one  after  deformation.  The  micro- 
densitometer is  used  to  digitize  the  film  negatives  by  assigning 
numbers  ranging  from  1 to  400  (-1%)  corresponding  to  light  intensity 
transmitted  through  the  negative  [49].  The  instrument  records 
approximately  300  density  units  per  second,  with  software  permitting 
incremental  step  sizes  between  1 pm  and  255  pm  (+2  pm)  [49].  A 
pattern  recognition  technique  is  then  used  to  obtain  displacements. 

The  method  ultimately  chosen  was  that  of  speckle  diffraction 
interferometry,  which  will  now  be  described. 


Speckle  Diffraction  Interferometry 

A brief  review  of  some  basic  concepts  of  light  will  define  the 
terms  needed  to  describe  the  interferometric  process.  Light  may  be 
regarded  as  an  electromagnetic  disturbance,  as  characterized  through 
Maxwell's  equations  [50].  A disturbance  moving  in  the  z -direction 
with  velocity  c has  the  form  [51]: 

A sin  (~  (z  - ct)  + 4>) 

p 

where  A is  the  amplitude  (A  is  called  the  intensity),  X the  wave 
length,  t the  time,  and  <j>  the  phase  angle  of  the  disturbance.  The 
receiver  of  the  electromagnetic  disturbance  (retina,  photographic 
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plate,  photomultiplier,  etc.)  responds  only  to  intensity  and  is 
insensitive  to  variations  of  phase  [521.  A continuous  series  of 
sinusoidal  cycles  is  called  a wave  train.  Two  wave  trains  are  said 
to  add  constructively  when  the  peaks  and  valleys  of  the  two  trains 
coincide  in  the  spatial  dimension,  and  are  said  to  add  destructi vely 
when  the  peaks  of  one  coincide  with  the  valleys  of  the  other. 
Constructive  addition  of  wave  trains  reinforces  and  increases  bright- 
ness (intensity), while  destructive  interference  tends  to  cancel  the 
disturbances  with  an  apparent  reduction  in  brightness  (intensity). 

When  a wave  train  is  divided  by  a beam  splitter  and  reunited,  it 
can  form  interference  fringes  provided  the  difference  between  the  two 
path  lengths  is  less  than  the  effective  length  of  the  characteristic 
wave  train.  Light  which  is  capable  of  interference  is  called 
"coherent",  and  the  characteristic  wave  train  length  is  called  the 
"coherent  length". 

The  coherent  length  of  lasers  is  several  feet,  while  that  of  a 
typical  light  source  is  on  the  order  of  1 mm  [51].  This  explains  why 
holography  did  not  become  an  important  tool  until  after  lasers  became 
readily  available  in  the  1960's.  An  excellent  review  of  holography 
is  given  by  Gabor  (531, while  a good  physical  introduction  to  holo- 
graphy is  given  by  Brown  [641.  Reviews  of  interferometry  can  be 
found  in  [55-57]. 

Consider  a diffuse  surface  illuminated  by  highly  coherent  light. 
Any  point  in  front  of  the  surface  receives  contributions  of  light 
reflected  from  all  illuminated  points  on  the  surface.  When  such  a 
surface  is  viewed  by  an  observer,  it  appears  to  have  a grainy  or 
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speckled  appearance  due  to  random  interference.  If  the  surface  is 
then  allowed  to  undergo  a displacement,  a new  speckle  pattern  is 
formed  corresponding  to  the  displaced  material  points  of  the  surface. 
If  a double  exposure  of  a photonegative  is  made,  and  the  negative 
then  developed,  we  have  a recording  of  the  two  speckle  patterns 
resulting  in  the  speckle  pairing  shown  in  Figure  15.  The  speckle 
size  is  proportional  to  the  wave  length  of  the  illuminating  light  and 
the  numerical  value  of  the  recording  device's  lens  aperture  [58]. 

The  governing  relationship  of  the  speckle  diameters  is  given  by  [59]: 


d = X F (1  + M) 


(26) 


where  x is  the  illuminating  light  wave  length,  F the  camera  f-  stop 
number,  and  M the  camera  lens  magnification.  For  the  helium-neon 
laser  used  in  this  work  (x  = 6328  X)  the  resolving  power  of  the 
emulsion  on  the  photographic  plate  had  to  be  at  least  1300/F  lines 
per  mm  [58,  59].  In  photographic  film  (plate),  the  higher  the 
resolving  power  of  the  film  (plate),  the  longer  the  exposure  time 

required  [60]. 

The  interrogation  of  a double  exposure  "specklegram" 
(photonegative)  consists  of  illuminating  the  specklegram  with 
coherent  light  from  the  laser  and  observing  the  light  intensity 
distribution  on  a screen  placed  so  as  to  collect  light  transmitted 
through  the  specklegram.  The  resulting  intensity  distribution 
collected  on  the  screen  can  be  shown  to  be  numerically  equal  to  the 
diffraction  halo  from  a single  exposure  transparency,  modulated  by 
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Young's  fringes  [61].  A schematic  of  the  interferometric  technique 
is  shown  in  Figure  16,  while  the  equipment  used  in  the  specklegram 
recording  and  interrogation  are  shown  in  Figures  17  and  18  respec- 
tively. The  specklegram  was  made  at  Texas  A&M,  while  the  interroga- 
tion of  the  specklegram  was  done  on  equipment  at  the  Flight  Dynamics 
Laboratory  in  Ohio.  Notice  that  recording  is  performed  on  an  optical 
steady  table  to  isolate  the  test  system  from  surrounding  noises. 

The  modulated  diffraction  halo  (Figure  19)  is  predicted  by  the 
Fraunhofer  diffraction  integral  which  is,  apart  from  a multiplicative 
factor,  a Fourier  transform  1611.  The  Fraunhofer  diffraction  integral 
applies  to  parallel  light  rays  and  is  a specialization  of  the  more 
general  Fresnal -Ki rchoff  integral  [55]. 

The  magnitude  of  the  displacement  field  of  the  specklegram  at 
the  point  of  interrogation  is  then  obtained  through  the  simple 
algebraic  equation 


6=  (n  -1/2 )x  M 

sin  « 


(27) 


where  n is  the  number  of  intensity  minimums,  M the  magnification  of 
the  specklegram,  \ the  wave  length  of  the  interrogating  beam  and  « 
the  fringe  minimum  angular  displacement. 

The  fringe  pattern  was  recorded  by  a linear  silicon  diode  (9 
inches  by  1 inch)  and  plotted  with  the  Varian  F80A  x-y  plotter  shown 
in  Figure  20.  Two  stepping  motors  moved  the  specklegram  past  a rigid 
laser  beam  in  a scanning  process.  Data  point  locations  on  the 
specklegram  were  recorded  with  an  accuracy  of  + .01  inches  (.02  inches 
on  the  specklegram  since  M =1/2). 
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A Mamiya  RB67  single  lens  reflex  camera  -'ith  a 90  mm  lens  was 
used  in  specklegran  i wordings . S pec k 1 eg ra were  recorded  on  Agfa 
Gavaert  10E75,  a high  resolution  ei.  on  on  a glass  plate.  The 
emulsion  thickness  on  the  plate  was  10  pm  with  a resolving  power  of 
2500  lines/mm.  It  can  be  shown  (Archbold,  Burch,  and  Ennos  [”58]) 
that  with  an  emulsion  thickness  of  10  pm,  an  f-stop  number  of  at 
least  2 must  be  used  so  that  the  speckle  pattern  does  not  change 
appreciably  with  depth;  an  f-stop  number  of  3.8  was  used  in  this 
work.  It  is  easily  shown  through  equation  (26)  that  a resolution 
of  277  lines/mm  is  the  minimum  resolution  tolerable  for  an  f-stop 
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number  of  3.8  and  light  wave  length  of  6328A. 

A 15  milliwatt  Spectra-Physics  124A  Helium-Neon  laser  was  used 
to  illuminate  the  specimen  with  a 10X  Bausch  and  Lomb  microscope 
objective  acting  as  a beam  diffuser.  The  relatively  low  power  of 
the  laser  restricted  the  area  of  illumination  on  the  specimen  to  a 
circle  with  approximately  a 2.5  inch  diameter.  After  double  exposure 
of  the  photographic  plate,  the  negative  was  developed  in  Kodak  D-19. 

Results 

Six  specimens  were  tested:  one  natural  unfilled  rubber,  one 
uncracked  propellant,  and  four  cracked  specimens  (two  .2  inches  thick 
and  two  .125  inches  thick).  A table  model  Instron  was  used  to  apply 
a constant  overall  strain  increment  with  nominal  grip  displaceme1 
of  .004  inches  for  each  increment.  The  total  overall  gr  p displace- 
ment of  .012  inches  was  reached.  A specklegram  was  made  corre- 
sponding to  each  applied  strain  increment  and  designated  according  to 
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specimen  and  grip  displacement;  e.g.,  6-2  was  the  designation 
accorded  sample  6 undergoing  the  second  strain  increment.  It  should 
be  noted  that  in  all  cases  the  first  applied  grip  displacement  gave 
ambiguous  readings,  perhaps  due  to  a healing  effect  between  the  two 
precut  crack  surfaces. 

The  rubber  sample  and  uncracked  propellant  samples  were  used  to 
determine  experimental  procedures.  The  rubber  sample  produced 
fringe  patterns  of  high  contrast.  Although  the  uncracked  propellant 
sample  produced  fringes  that  were  detectable  to  the  human  eye;  there 
was  not  sufficient  contrast  for  use  with  the  semi -automated  interro- 
gation equipment. 

An  examination  of  the  stressed  unpainted  propellant  sample  of 
Figure  1 (p.  2)  indicates  a large  amount  of  microcracking,  which 
causes  decorrelation  due  to  relatively  large  local  displacements  and 
particle  rotations. 

As  previously  mentioned,  the  interrogating  process  is  accom- 
plished by  passing  a lmin  laser  beam  through  the  speck! egram.  The 
fringe  pattern  produced  is  due  to  all  speckle  pairs  illuminated  by 
the  interrogating  beam.  The  fewer  the  number  of  correlated  pairs, 
the  weaker  the  fringe  pattern  contrast. 

PBAN  propellant  was  found  to  exhibit  significant  decorrelation 
outside  the  displacement  range  of  1-100  pm.  When  a natural  unfilled 
rubber  was  tested,  the  upper  displacement  limit  was  extended  to 
250  pm.  At  the  Flight  Dynamics  Laboratory,  under  ideal  conditions,  a 
pure  translation  of  750  urn  has  been  recorded. 

As  previously  mentioned,  the  effects  of  microcracking  and 
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particle  rotation  were  overcome  by  the  application  of  a thin  film  of 
flat  white  Krylon  spray  paint  applied  to  the  sample  surface.  The 
effect  of  the  paint  film  on  the  displacement  field  was  investigated. 
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To  test  the  effect  of  the  paint  on  the  displacement  field,  the 
upper  half  of  the  uncracked  specimen  was  masked  before  painting.  After 
painting,  the  mask  was  removed,  leaving  a distinct  paint  line.  A 
specklegram  of  the  specimen  under  load  was  made.  A scan  across  the 
paint  line  in  the  specklegram  proved  the  paint  film  had  negligible 
effects.  This  was  to  be  expected,  since  the  interrogating  laser  beam 
has  an  averaging  effect  over  the  area  of  the  beam  (1mm  diameter)  and 
the  effect  of  the  paint  was  to  average  displacements.  The  averaging 
effect  of  the  paint  can  be  seen  in  the  SEM  study  shown  in  Figures 
21  and  22. 


Discussion  of  Resul ts 

A comparison  between  the  numerically  obtained  displacement 
field  and  the  measured  displacement  field  was  made.  A comparison 
of  strain  was  attempted,  requiring  a differentiation  of  the  dis- 
placement fields.  A two-dimensional  cubic  spline  technique  was  used 
in  an  attempt  to  obtain  the  strains  associated  with  the  measured 
displacement  fields.  Data  scatter,  albeit  nominally  within  6%  of 
theoretical  results,  was  apparently  not  small  enough  for  the  success- 
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ful  use  of  the  spline  technique.  It  was  felt  data  scatter  was  due 
in  part  to  the  facilities  used  to  run  the  experiment  and  in  part  due 
to  local  failure  in  the  propellant,  i.e.,  large  movement  in  the  paint 
caused  by  paint  fracture  (the  interferometric  technique  has  dis- 
placement measurement  accuracies  of  1 micron  - 2%  [62]).  In  order  to 
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check  on  the  accuracy  of  the  spline  technique  itself,  in  fracture 
problems,  strain  fields  associated  with  rigid  body  displacements, 
affine  displacements,  Williams'  singular  crack  problem,  Barenblatt's 
singular  crack  problem,  and  the  combined  Wi 1 1 iams-Barenblatt  problem 
were  calculated  from  theoretical  displacements.  All  were  in  agree- 
ment with  the  analytically  derived  strain  fields. 

Of  the  four  cracked  specimens,  only  specimen  6 (.2  inches  thick) 
and  specimen  8 (.125  inches  thick)  had  fringe  patterns  with  suf- 
ficient contrast  to  interrogate.  The  other  two  specklegrams  (each 
cracked  sample  was  tested  only  once  due  to  material  damage  suffered 
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by  propellant  under  strain)  had  poor  contrast  due  to  overexposure 
and/or  over  development. 

The  specklegram  interrogation  revealed  still  another  problem. 

The  crack  tip  neighborhood  enclosed  by  the  lines  x=-0.1,  x=0.2, 
y=0.0,  y=0.1  was  consistently  absent  of  fringe  patterns  in  all 
specklegrams. 

There  are  seemingly  two  features  in  the  interferometric 
technique  which  may  be  violated  in  this  neighborhood,  causing 
speckle  pair  decorrelation:  time  dependent  motion  of  the  specimen 
and/or  rotational  movement  of  the  specimen  (paint  cracking  was  not 
considered  a problem  since  paint  flakes  were  considered  large 
enough  to  interrogate,  e.g.,  those  shown  in  Figures  21  (p.  54)  and 
22  (p.  55)). 

It  can  be  shown,  through  the  governing  field  equations,  that 
for  an  incompressible  elastic  body  with  displacement  boundary 
conditions  on  part  or  all  of  the  boundary  and  zero  tractions 
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specified  over  the  remaining  portion  of  the  boundary,  that  the 
strain  and  displacement  distributions  are  independent  of  the  modulus; 
therefore,  in  the  nonaging  linear  viscoelastic  case,  the  cor- 
respondence principle  implies  that  a similarly  loaded  viscoelastic 
material  has  strain  and  displacement  distributions  which  are 
independent  of  the  relaxation  modulus.  As  previously  indicated, 
that  portion  of  the  quarter  biaxial  strip  boundary  located  between 
the  apparent  crack  tip  and  the  crack  tip  is  acted  on  by  a non-zero 
stress,  which  implies  the  displacements  are  modulus  dependent. 

Rather  than  directly  treating  the  problem  of  the  non-zero 
boundary  tractions,  the  second  limiting  feature  of  the  interfero- 
metric technique  was  investigated. 

While  pure  translatory  motion  can  be  measured  over  a consider- 
able range  by  the  interferometric  technique,  the  degree  of  rotational 
movement  that  can  be  tolerated  is  more  severely  limited  [63]. 

Archbold  and  Ennos  [63]  record  fringe  patterns  with  rotations  of 
.02  radians. 

A scanning  electron  microscope  (SEM)  study  was  performed  on  a 
cracked  specimen  under  5%  loadings  to  qualitatively  determine 
paint  flake  rotations.  These  strains,  5X,  though  considerably 
higher  than  those  applied  to  the  test  specimen  (.2%),  were  the 
lowest  that  could  be  applied  using  the  relatively  insensitive 
tensile  stage  (Figure  23)  furnished  with  the  JOEL  JSM  - U3 
scanning  electron  microscope.  The  samples  were  ion  sputter  coated 
with  60:40  Gold  Paladium  alloy  to  insure  good  secondary  electron 
imaging  and  to  minimize  charging  effects  of  the  electron  beam  on 
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j the  sample. 

The  study,  as  recorded  in  Figures  21  (p.  54)  and  22  (p.  55), 
indicates  paint  flake  rotations  are  appreciably  reduced  away  front 
the  crack  tip  after  a distance  of  .01  inches.  Since  the  SEM  sample 
was  strained  approximately  25  times  higher  than  the  interferometric 
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COMPARISON  OF  THEORY  AND  EXPERIMENT 

All  grip  displacement  increments  applied  to  the  samples  were  nomi- 
nally .004  inches.  To  better  visualize  experimental  grid  deformations . 
these  results  were  normalized  to  a grip  displacement  of  .1  inches. 

Figure  24  compares  the  singular  crack  solution  as  obtained  through 
TEXGAP  (where  the  failure  zone  length,  a,  is  equal  to  zero)  with  nor- 
malized experimental ly  obtained  data.  Figures  25  and  26  compare  the 
experimentally  obtained  data  to  finite  stress  solutions  having  failure 
zone  lengths  of  .1  inches.  The  elements  generated  by  the  TEXGAP  code 
are  shown  in  Figures  24-26  as  solid  vertical  and  horizontal  lines. 

The  theoretical  solutions  shown  in  Figure  25  were  obtained  by 
superposing  previously  derived  Barenblatt  solutions  for  the  cracked 
infinite  sheet  problem  for  both  the  constant  and  ramp  failure  stress 
distributions  to  Williams'  singular  solution  (TFXGAP  solution)  for  the 
quarter  biaxial  strip. 

Figure  26  illustrates  the  TEXGAP  solution  obtained  by  linking  the 
failure  stress  distribution  with  the  singular  Williams'  solution 
through  the  finite  stress  condition  (Appendix  B). 

The  trend  throughout  most  of  the  strip  is  that  of  higher  experi- 
mentally measured  vertical  displacements  relative  to  theoretically 
predicted  displacements.  This  discrepancy  is  more  readily  recognized 
in  Figures  27  through  32,  which  expands  the  vertical  scale  by  a factor 
of  10  and  displays  the  displacements  of  Figure  26  (with  the  addition  of 
test  run  8-2).  There  were  no  horizontal  displacements  recorded  because 
a horizontal  shift  of  the  loading  frame  when  loads  were  applied  com- 
pletely masked  specimen  horizontal  displacements.  The  horizontal  load- 
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Figure  24.  Williams'  grid  deformations  (solid  curves)  for  the  quarter 
biaxial  strip  compared  to  experimentally  obtained  results 
( o contains  data  points). 
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Figure  25.  Williams'  grid  deformations  (TEXGAP  solution)  for  the  quarter 
biaxial  strip  superposed  with  Barenblatt's  grid  deformations 
for  the  infinite  sheet  compared  with  measured  deformations 
( o contains  data  points)  for  both  constant  (dashed  curves) 
and  ramp  (solid  curves)  failure  stresses. 


v (inches) 


Modified  Williams'  displacements  for  the  quarter  biaxial  strip  compared 
with  measured  displacements  0 y = .1  inches  with  the  vertical  scale 
expanded  by  a factor  of  10. 
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Figure  28.  Modified  Williams'  displacements  for  the  quarter  biaxial  strip 
compared  with  measured  displacements  (a  y = .2  inches  with  the 
vertical  scale  expanded  by  a factor  of  10. 
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Figure  30.  Modified  Williams'  displacements  for  the  quarter  biaxial  strip 
compared  with  measured  displacements  P y = .6  inches  with  the 
vertical  scale  expanded  by  a factor  of  10. 
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ing  frame  shifts  were  between  5%  and  10%  of  applied  grip  displacements. 

Four  possible  causes  of  the  relatively  high  vertical  displace- 
ments were  investigated:  1)  sample  thickness  effect,  2)  different 
Poisson's  ratios,  3)  different  failure  zone  lengths,  and  4)  nonlinear 
material  behavior. 

From  a strictly  theoretical  point  of  view,  one  would  suspect,  all 
else  being  equal;  the  thinner  the  biaxial  strip,  the  more  closely 
approximated  the  plane  stress  hypothesis.  This  contradicts  the  experi- 
mental results  recorded  in  Figures  27-32  (pp.  64-69).  Notice  in 
particular  for  y = .1,  .2  inches,  the  trend  of  the  thinner  sample's 
(test  run  8-2)  displacements  is  definitely  higher  than  the  displace- 
ments recorded  from  the  thicker  samples  (test  runs  6-2,  6-3);  however, 
as  we  move  farther  away  from  the  crack  centerline,  we  record  a closer 
agreement  between  the  thinner  .125  inch  thick  material  (test  run  8-2) 
and  the  theoretically  obtained  curve  than  between  the  thicker  .2  inch 
thick  material  (test  runs  6-2,  6-3)  and  the  theoretically  obtained 
curve.  Thus,  the  recorded  phenomenon  cannot  be  explained  as  a sample 
thickness  effect. 

It  is  well  established  that  Poisson's  ratio  for  propellant  is 
approximately  equal  to  .5,  i.e.,  propellant  behaves  as  an  incompress- 
ible material.  Propellant  castings  can  contain  large  internal  voids, 
thereby  exhibiting  compressible  behavior  with  effective  Poisson's 
ratios  substantially  lower  than  .5.  Accordingly,  Poisson's  ratios  of 
.45,  .4,  .3,  and  .2  were  used  in  additional  calculations.  The 
resulting  theoretical  displacements  in  both  the  horizontal  and  verti- 
cal directions  were  nominally  10%  smaller  for  a Poisson's  ratio  of 
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.45,  with  a similar  trend  continuing  (an  additional  4%  displacement 
reduction  for  each  reduction  in  Poisson's  ratio)  as  smaller  Poisson's 
ratios  were  tested.  We  can  conclude  that  the  high  vertical  displace- 
ments were  not  due  to  a Poisson's  ratio  effect. 

Finally,  an  investigation  of  the  effect  of  varying  failure  zone 
lengths  was  conducted.  Figures  24  (p.  61)  and  26  (p.  63)  illustrate 
the  cracked  biaxial  strip  with  failure  zone  lengths  of  0 and  .1  inches, 
respectively.  The  failure  zone  length  was  felt  to  be  less  than  .2 
inches  because  of  the  speckle  decorrelation  effect  mentioned  earlier. 

A failure  zone  length  of  .2  inches  (not  shown)  was  tested.  The  trend 
established  by  Figures  24  (p.  61)  and  26  (p.  63)  continued;  i.e.,  the 
effect  was  limited  to  the  crack  tip  neighborhood  with  the  theoretically 
predicted  vertical  values  higher  than  those  shown  in  Figure  26.  For 
all  of  the  above  failure  zone  lengths  the  measured  vertical  displace- 
ment field  away  from  the  crack  tip  remained  high.  Evidently,  the 
high  experimentally  obtained  vertical  displacements  cannot  be 
theoretically  explained  by  varying  failure  zone  lengths. 

The  sensitivity  of  the  experimental  technique,  as  stated  pre- 
viously, is  - 1 micron,  which,  for  the  displacements  used  in  this 
experiment,  represent  an  accuracy  of  approximately  2%  — not  suffi- 
ciently large  enough  to  account  for  the  observed  phenomenon. 

We  are  left  to  conclude  that  the  displacement  field  behavior 
cannot  be  explained  through  linear  theory,  and  must  represent  a 
nonlinear  behavior. 


i 


72 


A QUESTION  OF  UNIQUENESS 
Motivation 

There  remain  several  open  questions  concerning  crack  problems, 
and  among  these  are:  1)  is  the  Barenblatt  solution  unique,  and  2)  are 
the  displacements  the  smoothest  we  can  hope  for,  i.e.,  ueC(fi)  with 
derivatives  failing  to  hold  throughout  the  problem's  domain,  n? 

Knowles  and  Pucik  [64]  address  the  question  of  uniqueness  of 
solution  to  the  traction  boundary  value  problem  for  the  two- 
dimensional  linear  elastostatic  problem  for  a bounded  domain  containing 
a crack.  In  addition  to  the  standard  hypotheses  in  the  Kirchoff 
uniqueness  theorem  (cf.  Love  [65]),  Knowles  and  Pucik  require  the 
displacement  field  to  be  bounded.  With  these  hypotheses,  Knowles  and 
Pucik  show  that  the  strain  energy  in  a neighborhood  of  the  crack  tip 
remains  finite,  and  then  prove  uniqueness. 

In  the  following  section  we  address  the  question  of  uniqueness 
of  solution  to  the  traction  boundary  value  problem  for  the  three- 
dimensional  linear  elastostatic  problem  for  a class  of  bounded  regions 
(that  will  be  made  precise  later  in  the  chapter)  by  relaxing  the 
Kirchoff  stress  and  displacement  smoothness  hypotheses.  As  a special 
result,  uniqueness  for  the  three-dimensional  crack  problem  is  obtained. 

The  uniqueness  theorem  of  the  next  section  does  not  preclude 
solutions  with  infinite  strain  energy.  The  stresses  need  only  be 
absolutely  integrable  in  the  Lebesgue  sense,  whereas  finite  strain 
energy  implies  square  integrabi 1 ity. 

Our  formulation  requires  that  ueC(n),  which  implies  bounded 
displacements  for  bounded  domains. 
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Classical  Elastostatic  Results 

We  begin  with  a review  of  some  classical  elastostatics  results; 
however,  before  we  begin  with  the  review,  some  necessary  notation 
will  be  described.  Ordinary  lower  case  letters  represent  scalars, 
while  letters  with  a tilde  placed  under  them  represent  tensors  of 
positive  order.  We  designate  e.y  as  the  fully  contracted  outer 
product  and  write 


^ = ciJYiJ 


where  repeated  indices  imply  summation.  Furthermore,  we 

define  ||  a>||  of  any  order  tensor  to  be  the  square  root  of  the  fully 

n 

contracted  outer  product  of  w.  We  write  R = II  R-  where  R is  the 

n 1=1  1 

real  line,  and  (Rn,U)  the  usual  topological  space,  where  U is  the 

topology  generated  by  the  projection  mappings.  By  a domain  in  Rn  we 

mean  an  open  connected  subset  in  Rn,  by  a region  we  mean  a connected 

subset  in  R . Let  ilcR  then  aft,  ft,  ft  denote  the  boundary,  interior 
n n 

and  the  closure  of  ft  respectively. 

Let  ftcRn,  then  C(ft)  is  the  set  of  all  tensor  valued  functions 
continuous  on  ft,  and  Cm(ft)  is  the  set  of  all  tensor  valued  functions 
continuous  on  ft  and  with  all  partial  derivatives  up  to  order  m,  a 
positive  integer,  existing  and  coinciding  with  functions  C(ft).  We 
further  write 


C (ft)  » O C"(ft). 

n»l 


1 1 d r field,  is  a function 
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|l 


L 


defined  on  n with  range  in  R.  We  write  ■*  R.  We  use  the  standard 
covariant  differentiation  symbols,  i.e.. 


m 


9 <$> 


’ i j . . . k 


ax... .ax. ax. 

k j i 


m indices 


Here  the  covariant  differentiation  raises  <j>  to  an  mth  order  tensor. 
Similarly,  let  u be  a first  order  tensor  in  C1  (S2) » then  we  write  vu 


for  u.  . where  u.  . is  a second  order  tensor. 

1 1 j J 


We  now  state  the  Divergence  Theorem  [66]  for  regular  regions, 
where  regular  region  is  an  adoption  of  Kellogg's  [67]  notation. 

Divergence  Theorem  - Let  n be  a bounded  or  unbounded  regular 
region.  Let  4>  be  a scalar  field,  u a vector  field,  and  T a tensor 
field,  and  let  <j>,  u,  and  T be  continuous  on  n,  differentiable  almost 
everywhere  on  n,  and  of  bounded  support.  Then 


'an 


'an 


an 


an 


<j>n.dS  = / 

1 

* ,d  V 
♦ 1 

u(„.dS  -J 

u,  .dV 

U.HjdS  -jf 

u . sd V 
i ,J 

cijk"jukds 

slc' 

T . .n .dS  = 
1J  J J 

^ Ij.J 

.d  V 


dv 


whenever  the  integrand  on  the  right  is  piecewise  continuous  on  n. 

We  now  state  some  classical  results  from  the  linearized  theory  of 
elasticity  (cf.  [68]). 


75 


Definition  (State:  reyular  elastostatic  state.)  Let  R be  a 
bounded  domain  in  R^,  and  DcR.  If  u is  a vector- valued  and  o a 
second-ordered  tensor-valued  function  on  D,  the  ordered  pair 
S = (u,g]  is  called  a state  on  D with  the  displacement  field  u and 
the  stress  field  o. 

Suppose  thi  t f is  a vector  field  on  D and  c a fourth  order 
tensor  field  on  R.  Then,  S = [u,g]  is  said  to  be  a regular  elasto- 
static state  on  D corresponing  to  the  body-force  field  f and  the 
elasticity  field  c,  provided: 

i)  frC(D) ; CeC(R)  with 

cuki  ■ cjiki  ■ ckiij on  R 

and  c invertible  on  R,  i.e.,  the  six  by  six  matrix  associated 
with  c has  nonvanishing  determinant  on  R. 
if)  grCZ(n)Oc'  (0) . orc'lo)  n C([>) 

if1>  “ij- s-jkl  uk,l  in  “ 

o...+f.  =0  in  D 
1J.J  1 

We  write  E(f,c;D)  for  the  set  of  all  regular  elastostatic  states  with 
the  above  f,  c,  and  D. 

Let  S'  = [u'.g'l  and  S = [u,g]  be  in  E(f,c;D)  and  let  cuR. 

We  define  addition  and  scalar  multiplication  of  two  states  in  E(f,g;D) 
through: 


and 


S + S'  = [ U ,g ] + ( u ' , o * | = [u  + u',  g + g'] 
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1 I.  mi  — -T— * 


L 


«S  = a[u,o)  = [aU,ao] 

implying  E(f,c;D)  is  a linear  vector  space.  Regular  regions  are 

composed  of  a finite  number  of  regular  surfaces,  which  may  have  in 

common  a single  point  (a  vertex  for  both),  or  a single  regular  arc  (an 

edge  for  both),  but  no  other  points.  At  these  common  points  there 

may  not  be  a definable  normal.  With  the  preceding  as  motivation,  we 

★ 

define  the  regular  part  of  DR,  designated  DR,  as  that  subset  of  DR 

at  which  a normal  is  definable.  Let  R be  a regular  region,  S = [u,o] 

a state  on  R,  and  n the  unit  outward  normal  on  dR,  then  the  traction 
★ 

vector,  s,  of  S on  dR  is  defined  by 

it 

s . = o . .n  • on  DR 
^ U J 

We  now  state  the  traction  problem  for  regular  elastostatic 
states  (classical  traction  problem).  Let  ncR  be  a regular  region. 

We  wish  to  find 


S = [u  ,o ] e E(f , c ; ^ ) 

given  the  data  f and  c subjected  to  the  boundary  conditions 

★ ★ 

s 3 § on  DR 

★ 

where  s is  a prescribed  function. 

The  uniqueness  of  the  solution  to  the  classical  traction 
problem,  due  to  Kirchhoff,  is  now  given  (cf.  Love  [65)), 

Theorem  (Uniqueness  theorem  of  classical  elastostatics: 
Kirchhoff).  Let  S2  be  a regular  region,  and  suppose 
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a)  S'  = [u'  ,o‘ ] , S"  = [u",a“]  c E(f,c;S) 

* * ★ 

a)  5*  = §"  on  3SJ 

y)  c is  positive  definite  on  n,  i.e., 

c1jkl ^Yij^Ykl ^ > 0 on  for  311 

and  for  every  non-zero  symmetric  two-tensor,  y.  Then  S'  = S"  on  R 
modulo  a rigid  displacement. 

Now  we  relax  the  previous  regularity  assumptions  on  the  solution 
to  the  traction  problem  of  elastostatics  and  deal  with  the  attending 
problem  of  uniqueness.  After  dealing  with  necessary  preliminaries, 
we  cast  the  traction  problem  in  its  proper  setting,  and  then  prove  an 
appropriate  uniqueness  theorem.  The  resulting  arguments  closely 
follow  the  spirit  of  Brockway  1681,  who  looks  at  the  uniqueness  of 
solutions  to  the  traction  problem  for  "weak  elastodynamic  states". 

Mathematical  Preliminaries 

The  following  requires  liberal  usage  of  functional  analysis 
results.  Functional  analysis  rests  heavily  on  two  pillars  of 
mathematics.  Topology  and  Algebra.  As  such,  the  following  definitions 
and  results,  albeit  not  complete,  will  aid  the  engineering 
reader  through  the  maze  of  terminology  [69-74]. 

A collection,  t,  of  subsets  of  a set  X is  said  to  be  a 
topology  in  X^  provided: 
i)  <f>et  and  Xcx 

ii)  If  O.ct  for  i=l,...,n;  then  O Oct 
1 i i=l  1 


■ 

I 


. TBH.- 
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iii)  If  0 ex,  for  all  ael,  then  u 0 ex. 
a ael  a 

Special  note  should  be  made  of  the  finite  number  of  intersections  in 

i i ) and  of  the  finite,  countable,  or  uncountable  number  of  sets  in 

iii) . 

A family  of  sets  is  countable  provided  there  exists  a bijection 

(one-to-one  correspondence)  between  a subset  of  the  positive  integers 

and  the  collection  of  sets,  i.e.,  we  can  label  the  sets  0^ , Ogi 

If  t is  a topology  in  X,  then  (X,t)  is  called  a topological 

space.  Let  (X,x  ),  (Y,x  ) be  topological  spaces,  and  f a mapping 
x y 

(function)  of  X into  Y . Then  f is  said  to  be  continuous  provided 
f' ~^(v)et  for  every  vex  . By  function  we  mean  "single-valued" 

X jr 

function  in  the  standard  mathematical  sense. 

Let  M be  a collection  of  subsets  of  the  set  X,  then  M is 
called  a o - algebra  provided: 

i)  XcM 

ii)  if  AeM  then  X\A  = {xrX|x{:A}eM 

on 

iii)  if  { A„}  is  a countable  number  of  sets,  then  u A eM. 

n n=l  n 

Let  M be  a o-algebra  in  X,  then  the  ordered  pair  (X,M)  is  called 

a measurable  space  and  the  members  of  M the  measurable  sets  in  X.  If 

(X,M)  is  a measurable  space  and  (Y,x)  a topological  space,  then  a 

mapping,  f,  of  X into  Y is  said  to  be  measurable  provided  f ^(v)i-M 

for  every  Vex. 

A positive  measure,  p,  is  a function  defined  on  a o-algebra, 

M,  whose  range  is  in  [0,«]  and  which  is  countably  additive,  i.e., 

1 f { A ^ } is  a countable  number  of  disjoint  sets  in  M (A^HAj  = <j> 

oo  oo 

whenever  i / j ) then  u(u  A.)  = Z p(A.).  To  avoid  trivial  cases  we 

i=l  i=l 
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will  further  assume  there  exists  at  least  one  A^M  such  that  p (A ) < 
We  define  a measure  space  as  the  ordered  triple  (X,M,u).  We  now 
state,  without  proof,  some  facts  of  positive  measure. 

Let  u be  a positive  measure  on  a o-algebra  M.  Then 

a)  mU)  = 0. 

b)  p(A^U...UAn)  = u(Aj)  + ...  + u(An)  if  A^ , ...,  An 

are  pairwise  disjoint  members  of  M. 

c)  ACB  implies  u(A)  _<  w(B),  A,BcM. 

oo 

d)  u(A  1 ->  y (A)  as  n •*»  where  A = u A , 

n n=l  n 

A^eM  and  A^cA^  whenever  i < j. 

oo 

e)  w(An)  u(A)  as  n ->  » if  A = HAn,  A^M 

n=l 

and  A.3A.  whenever  i < j 
■ J 

and  u(Aj ) < ® . 

Let  (X,M)  be  a measurable  space.  The  function  s : X -*■  [0,®)  such  that 
the  range  of  s consists  of  only  a finite  number  of  points  in  [0,®) 
is  called  a simple  function.  Suppose  aj,...,an  are  the  distinct 
values  of  a simple  function,  s.  Set  Ai  = {x | s (x ) = a.},  define  XA 
through 


XA  (x)  = 1 x t A^ 

= 0 otherwise 


we  characterize  s through 

n 

S ( X ) = Z a.Xfl  (x) 
i=l  1 Ai 


and  denote  X.  as  the  characteristic  function  of  A.,  clearly,  s is 
Mi  1 


measurable  if  and  only  if  is  measurable. 

We  are  now  in  a position  to  define  Lebesgue  integral . 

Let  (X,M)  be  a measurable  space  and  s be  a measurable  simple 
function  on  X,  we  write 


n 

S = Z a.X 
i=l  1 1 


further,  let  EeM,  then  define 


/ 


n 


sdu  = 2 a-w(A.nE) 
E i=l  1 1 


We  use  the  convention  0‘<*  =0  (note,  some  might  be  0 and 
p(A.nE)  = » 

If  f : X ->  [0,«]  is  measurable,  and  EeM,  we  define 


/ fdu  = sup  / sdu 
JE  0<s<f  JE 


Occasionally  it  is  desirable  to  define  the  integral  for 

measurable  functions,  f,  with  range  in  [-<»,<») , the  obvious  result  is, 

(provided  at  least  one  of  the  integrals  / f+du»  J f du  is  finite) 

■'E  E 

fdM  = £f+d u 

where  f+(x)  = max  { f (x ) ,0 } and  f’(x)  * max  {-f(x),0}. 

An  easy  result  to  show,  for  non-negative  f and  s is 

<#•  ( E ) =/  sdu  and  $(E)  =/  fdu  define  positive  measures  on  M. 

J E JE 

Let  n be  a bounded  subset  of  R^,  and  0 < p < <*>.  Define  Lp(g)  by 


Lp(n)  = (f|f  is  tensor  valued  and  measurable  on  n and 
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Jr  |f..  . | Pdp  < oq  where  dp  is  the  Lebesgue  measure  in  R^}. 

^ 1 J . • . K o 

We  must  now  turn  our  attention  to  regular  surfaces  in  order  to 
clarify  the  divergence  theorem  and  extend  the  classical  (Riemann) 
integrals  to  Lebesgue  integrals. 

A surface  is  a mapping  a : n ■*  R3  where  ncR^.  The  trace  of 
a surface  is  the  range  of  8.  In  our  work  we  deal  with  regular 
surfaces  and  regular  regions  in  the  sense  of  Kellogg  (Kellogg's 

surfaces  coincide  with  the  above  definition  of  trace).  As  such, 
given  a regular  surface  S,  there  exists  a Cartesian  frame  in  R^ 

such  that  S admits  the  representation 

x3  = G(x1,x2)  for  all  (x^x^ell,  GecV) 

where  n is  a closed  region  in  the  (x.j  ,x2) -plane.  By  an  integral 
over  the  surface,  S,  of  a function,  we  mean 


where  f = (o,o,f)  and  n is  the  unit  normal  to  the  surface  and 
JG(xi  .x2)cC(nK  anc*  Jq(x)  is  the  Jacobian  of  G at  x (cf. 

Buck  [72]  , Rudin  [71  ] ) . 

We  need  more  mathematics,  consider  R^,  then  a set  of  the  form 

W = {xeRk|awCi-<3i , 1 < i 1 k where  x = ve{<,<_}} 

is  called  a k-cell.  Its  volume  is  defined  to  be 

k 

vol(W)  = n (p.  - a.) 
i=l  1 1 
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We  call  a set  of  the  form 

Q ( a *, 6 ) = {XeRk|a.  < Z-  < a i + 6 , 1 < 1 < k} 

a 5-box  with  a corner  at  a = (a.,..., a ). 

i n 

Let  (x,  M,  p)  be  a measure  space,  with  p a positive  measure. 

We  then  say  (x,  M,  p)  is  complete  or  p is  a positive  complete  measure 
on  M provided  given  any  AeM  such  that  p(A)  = 0,  and  EcA  then  EeM. 

We  now  state  a familiar  result  of  real  analysis  [71].  There 
exists  a positive  complete  measure,  p,  defined  on  a o-algebra,  M,  in 
Rk,  with  the  following  properties: 

i)  p(W)  = vol  (W)  for  all  k-cells  W 
ii)  p is  translate  invariant,  i.e. 

p(E  + x)  = p ( E ) for  all  EeM,  xeRk> 

We  call  the  members  of  M the  Lebesgue  measurable  sets  in  Rk  and 
p the  Lebesgue  measure.  We  now  state  the  divergence  theorem  in  a 
relaxed  setting,  the  proof  of  which  is  given  in  Appendix  C. 

Let  n be  a bounded  regular  region  (hence  a Lebesgue  measurable 
set),  and  $,  u,  T a scalar,  vector,  and  second  order  tensor  field, 
respectively,  differentiable  almost  everywhere  on  n with  derivatives 
in  L \n).  Then  the  results  of  the  classical  divergence  theorem  hold 
on  n with  the  surface  integral  interpreted  as  defined  on  page  81. 
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Uniqueness  of  Solution  to  the  Traction  Problem 
for  Weak  Elastostatic  States 

Let  n be  a bounded  domain  in  R . Define  Cn(S)  = 

n 

{feCn(«)|suppfcn},  where  suppf  = fxen|f(x)  f 0}.  Clearly,  if 

feC  (si)  then  f vanishes  in  a neighborhood  of  3fi. 

Definition  (generalized  gradient):  Let  f be  a vector  field, 
1,  . 

feL  (n)  where  n is  a bounded  domain  in  R^.  If  there  exists  a two- 

1 -l  _ 

tensor  field  weL  (n)  such  that  for  every  <(>cC  (fi) 
fa  Vij  i »je{l , 2,  3} 


we  say  f has  generalized  gradient  w on  n and  we  write 


3f  = u 

Suppose  there  exists  and  w"  in  L^fi)  that  satisfy  the  above 
definition.  Then 


toV.)  for  all  *cC1(n). 

I J 


Therefore  $•(<*>]•  - w"  ) = 0 almost  everywhere  (a.e.),  i.e.,  zero 

I I J I J 

except  on  a set  of  measure  zero.  The  arbitrariness  of  $ dictates 
. = uV.  almost  everywhere. 

* J I J 

Now  suppose  qeC  (6).  We  wish  to  show  vq  = 3g.  Applying  the 
divergence  theorem  to  the  first  integral  on  the  right  we  get 
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but  f <J>  .q  .n  . = 0 since  $eC  \n)  , 

“'an  1 J 

and  by  the  uniqueness  of  aq  the  result  follows. 

Definition  (W^si)):  Let  n be  a bounded  domain  in  R^, 

p _ 

X = {qeL  (n)|q  is  vector-valued,  aq  exists},  and 

«2<5)  ■ ( x.  II- II  s) 

where 

11911,  = /<l9l2  * l»9l2). 

n 

Remarks 

1.  X is  a vector  space  over  the  real  field,  where  the  elements 
of  the  space  are  sets  of  functions  that  are  equal  except  on  sets  of 
Lebesgue  measure  zero. 

2.  || • ||  is  a norm,  i .e. 

i)  ||  q||^  >_  0 and  ||  q||s  = 0 iff  q = 0 a.e.  for  all 

gcW^n) 

11)  II  q + p|ls  i II  q|ls  + II  Plls  for  all  p,  qew](S) 

iii)  II  oq  II  s = | a | ||  q ||  s for  all  qcWj(n),acR 

3.  w^n)  is  a complete  normed  linear  space,  i.e.,  every  Cauchy 
sequence  in  the  space  converges.  { f n > is  a Cauchy  sequence  provided 
given  any  e > 0 there  exists  MeP,  the  positive  integers,  such  that 
m,  n > M implies  ||  f - f ||  < e.  A complete  normed  linear  space  is 
called  a Banach  Space  (the  proof  of  completeness  is  established  by 

the  standard  argument  of  selecting  a candidate  in  for  a sequence's 


I 
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limit  and  then  showing  the  convergence  of  the  sequence  to  the 
limit,  cf  Smirnov  [75] 

To  motivate  the  definition  of  weak  elastostatic  state  about  to 
be  introduced,  consider  the  following.  Let  ft  be  a regular  region  and 

S = [u,o]  eE(f,c;5), 


We  wish  to  show 


or 


+j (a  v,  ■ 

/ (-V<}>.a  + <i>.f)  + / <t>.s  = 0 
« ' ' */3ft  ' ' 


for  all  <|>  e (ft) . 


(28) 


Let  p.(x)  = - # . (x)o • • (x)  for  all  xeft.  The  smoothness  of  p 

J * * J * 

impl ies 

But  SeE(f,c;ft)  (the  stress  equation  of  motion)  implies 

The  smoothness  of  p.  . allows  us  to  apply  the  classical  diver- 
J > J 

gence  theorem, 
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I 

I 

I 

I 

I 

* 

I 


~ — 


the  result  follows.  | 

Definition  (Weak  elastostatic  state).  Let  fi  be  a regular 
region,  and  c as  defined  before,  and  f a vector  field  measurable  on 
«.  Then  S = [u,g]  is  said  to  be  a weak  elastostatic  state  on  ft 
corresponding  to  the  body-force  field  f and  elasticity  field  c, 
provided: 

1)  ftL2(S),  ue  w’(5) 

i()  ■ cijkl<5)3luk(?) 

i i i ) f (-4>,  ;0..  + <t>.f.)  = 0 for  all  4>eC  (n) 
i.j  ij  ii 

For  the  set  of  all  such  states,  we  write 


Ew(f,c;6) 

What  is  the  relationship  between  E(f,c;6)  and  E (f ,c;fi)  ? 

' ■*  w ~ ~ 

Proposition:  E(f,c;fi)  cE..(f  ,c;fi) ; furthermore  if  SeE.(f,c;fi)  with 

' - W ~ ~ W ~ ~ 

p 1 1 

ueC  (n)HC  (n),  and  oeC  (ft)nC(n),  then  SeE(f,c;ft)  where 

f |f  - f'|2=  0,  i.e.  ||  f - f*  ||  2 = 0 
•A?  L c 

Proof:  E(f ,c;n)cE..(f ,c;fi)  follows  trivially  from  the  smoothness 

— - w **' 

condition  of  E(f,c;n)  and  equation  (28).  Now  let  ScE  (f,c;n)  together 
with  the  hypothesis  as  stated  above. 

Define  f1  through 

f‘  = - v.o  (f!  = - a..  ■ ) on  n 

1 1 J * J 

then  clearly 

S = [u,o]  EE(f ’.c’;  n)cEw(f ' ,c;ft) 
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f = f a.e.  iff  ||  f - f'|l2  - 0. 

\r 


Theorem:  Let  fi  be  a bounded,  locally  star-shaped  domain  in  R^.  Then 
C \n)  is  dense  in  W^(a). 

Proof : (Smirnov  [75]  proves  this  for  bounded  star-shaped  domains, 
this  proof  is  due  to  Brockway  [68]). 

NOTE:  If  ( X, t ) is  a topological  space,  then  A c X is  dense  in  X 
provided  A = X.  Here,  1^(6)  together  with  the  usual  topology  defined 
through  ||  • ||  is  a topological  space*  indeed  it  is  a Hilbert  space 
and  we  need  only  show  given  any  peW^)  that  there  exists  a sequence 
{pn)cc(n)  such  that  pn_!Llil^p  i.e.  ||  pn  - p||s  + 0 as  n + ® . 

If  xeQ,  then,  there  exists  a neighborhood  N(x)  such  that 

nnN(x)  is  star-shaped.  Now  (R^.ll  . ||  2)  1S  a finite  dimensional 

normed  linear  space  (note:  ||  x||  9 = x2  + x,  + x~). 

' ir  1 c J 
We  state  the  following  result  [73]. 

Proposition:  If  (V , || - 1|)  is  a normal  linear  space  over  the  field  $, 

then  the  following  are  equivalent: 

i)  (V, ||*||)  is  finite  dimensional 

ii)  If  EcV  is  a closed  bounded  set,  then  E is  compact. 

Definition:  Let  (X,t)  be  a topological  space,  then  A c X is 

compact  provided  every  open  cover  of  A is  reducible  to  a finite  cover. 


r! 


1. 


R8 


Let  r = {Gj  |G^et  ,icl) , we  say  r Is  an  open  cover  of  AcX  provided 
iVlGf^A*  ^e  Precoc|ing  proposition,  fi  is  compact,  therefore  there 
exists  a finite  subcovering  from  r = {N(x)|xrfi},  namely, 

{N(x’),N(x2) N ( xm ) } 

such  that 

m k - i 

UN(x  ) Dfi  and  such  that  N(x  ) is  not  contained  in 
k=l  ' 

m . 

U N(xJ). 
j“l  ' 

J*1 

Next,  let  = N(xk)  k = 1,  2,  ....  m. 

Define 


Dk  = Nkn«  k = 1 , 2,  . . . , m. 

Now  each  Dk  is  star-shaped.  If  peW^Ol)  then  pk  is  in  for 

all  k = 1 ...  .,m  where  pk  is  the  restriction  of  p to  Rk. 

Then  by  Smirnov  there  exists  {pk>cc  \f)k)  such  that 

P^LikkJ?  as  n -*•  » ( || • || sk  is  the  restriction  of  ||.||s  to  f>k). 

Define  H(fi)  = {feC“’(fi)|f  has  compact  support  in  fi}. 

We  state  the  following  well  known  theorem  [70]; 

Theorem:  If  A is  a collection  of  open  sets  in  Rn,  whose  union  is  fi, 

then  there  exists  a sequence  { . } c i:(fi)  with  ^ ' 0,  such  that 

«)  each  has  its  support  in  some  member  of  A, 

3)  2 = 1 for  all  xcfi 

i=l  1 

y)  to  every  compact  Kcfi  corresponds  an  integer,  m,  and  an  open 


I 


set  W c Ksuch  that 


'J’j(x)  + il^?)  + •••  + ^(x)  = 1 for  -tW" 

Such  a collection,  { <|> ^ > is  called  a locally  finite  partition  of  unity 
in  n subordinate  to  the  open  cover  A of  n. 


Corollary:  If  A is  a finite  collection  of  open  sets  in  the  above 
theorem,  then  there  exists  a sequence  with  ij.  >_0  such 


a)  each  ^ has  support  in  some  member  of  A 

b)  X ♦ * = 1 for  all  xifi,  where  there  are  n open  sets  in  A. 
i=l  1 

Proof:  Let  {Dj,...,Dm>  be  the  open  cover  of  n.  Then  corresponding  to 
each  there  exists  a compact  set  such  that  K.c  D..  Then  from  the 
above  theorem  there  exists 

\p-| i + ...  + . = 1 for  all  XeW.uKj,  W^..  open  and  W^Cfl. 

pi' 

At  least  one  of  the  if>..  has  support  in  D.  since  there  exists 

_ J 

m 

xcD.,x^u  D.  and  by  choosing  K.  such  that  xcK..  Now 
1 j=l  J 1 ' 1 

j=i 

i m Pi 

— X X 4>..  = 1 and  by  defining  ij».  = x<».  . such  that  supp  <K.c  D.. 

m -j  = ] j - ] *J  I IJ  J*  1 

The  result  follows.  ■ 


Introduce  a partition  of  unity  {i|»j  , . . . subordinate  to 
IDpD^,...  ,Dm),  where 

m 

E ^(x)  = 1 for  a1^  ancl  SUPP  l*'kcDk  ^ ~ 


90 


Define  q£  = i/^p"  on  (no  sum  on  k) 

= 0 on  (k=l . . ,m;n=l  ,2,. . . ), 

clearly  q£  = cC1 (fl). 

m 

Define  pn  = 2 q.  on  n , 

k=l  ~k 

clearly  p^C1^)  (n  = 1,  2,  ...). 

We  need  only  show 

n II  * II  s 

P — -P 

m m m n 

Now  II  P - P lls  =11  ^2  *k  - 9k H$  * \P  ' s’ 

but  ||  *kp  - qJJUs  = ||  H-k(p  - p")||Sk 

since  ^.g"  = 0 on  n\Dk  and  qj  = »PkpJJ  on  n. 

£ 

Now  ||  <Pk(p  - pj)||  Sk  -f  I |<Pk(p  - pj)|2  + I 3{^k(p  - p")  > 1 21 

Dk 

Simirnov  shows  the  product  rule  applies  to  yeneralized  differentiation 

1 CO  _ 

of  functions  in  W,,  . Since  i^eC  (fi),  i/>k  and  are  bounded  on  n, 
therefore  there  exists  MeR  such  that 

II  *k(p  - pjj)  ||  2 < M2/ {|p  - pj|2  + |3(p  - Pk|2l  for  all  k=l,...,m. 

k Dk 

Therefore  ||  i>k(p  - Pk)||s  ±M||  p - pj||  s 

k k 

But  II  P - pjll  sk  + 0 as  n +•  , 


which  implies  ||  p - pn||s  + 0 as  n -*• « . 

Therefore  prC^(ii)  and  C^(sl)  dW^(q)  . 

But  by  definition  of  W^fi),  W^nJoC^ii),  the  result  follows.  | 
We  formulate  the  weak  elastostatic  traction  problem.  Given 
a regular  locally  star-shaped  region,  n,  find 


s = f y , o ] r Ew  ( f ,c-,n) 


corresponding  to  data  f,  c and  subjected  to  the  generalized  boundary 


conditions 


I * ♦ffi)  '/  ■ 8 for  a"  »£C'(5> 


where  s is  L^(an)  as  discussed  previously  (p.  81). 

That  this  is  a generalization  of  the  classical  elastostatic 
problem  follows  by  assuming  enough  smoothness  on  S. 


Theorem:  (Uniqueness  of  solution  for  the  traction  problem  of  weak 
elastostatic  states):  Let  « be  a regular  locally  star-shaped  region 


а)  S ' = [u'.a’l  c Ew(f,c;a) 
S'1  - [Un  ,a"  ] e Ew(f,c;q) 

б)  for  all  $eC ^(si ) 


+ *fm  Vi 

ifi>  *L  V> 


= 0 where  §>  L^(as)) 


I 


Furthermore,  assume 

y)  c is  positive  definite  on  n,  so  that 
cijki(x)  Y^fxJYk^x)  > 0 for  all  xcfl 

for  all  symmetric  two-tensors,  y t 0.  Then  u'  = u"  modulo 
a rigid  displacement,  and  o'  = g"  almost  everywhere  on  n. 

Proof:  Define  a state  S =[u,g]  through  u = u'  - u",  g = o'  - g"  on  fi 

according  to  a),  b)  and  the  definition  of  weak  elastostatic  state, 

S = [y,g]eE  (0,c;fi),  f ,0..  = 0 for  all  feC1  (fi) . (A) 

w * » J 1 J 

We  show  that  (A)  holds  for  <t>ewl(n)  with  a.  . replaced  by  a . 

- C 1 »J  J • 

If  $ew!,(r:),then  since  C1  («)  is  dense  in  W^(n),  there  exists  a 
net  {$n}  c C1  (fi)  such  that  $nJil!k_$  a.e.  on  n. 

Define  a through 

a -fa^ i - • 

Now 

w - ♦? j>i  -jfhji  vi  -*"ji  • 

and  an  application  of  Holder's  inequality  implies 

|A|  - {jj ( l°ij|2}V2  {/  |aj*i  ' ♦i,j,2}1/- 

But  $,  $n  F-W^(il) , and  ||  An  - a||s  ■*  0 as  n -+•*  , i.e. 

- 4>n | 2 + 1 3$  - v^n | 2 } + 0 as  n + « . 

* 

Therefore  f | a>J>  - v^>n  | ^ + 0 as  n + » 

Jo 


and  | A | + 0 as  n >®.  That  is  to  say,  given  there  exists 


^ncC1(fi)  such  that 


But  (A)  implies  f .o..  - 0 for  all  Ac1  (fi)  , 
✓ ^ * J 

f * 0 . 

J = ° for  a11  , 

in  particular,  UeW^(n)t  therefore 

fau  jui)oij = °- 

But  S = (u,o]  is  a weak  elastostatic  state,  therefore 


therefore 

Now 


o..  = C - ■ . , 3,  u. 

ij  l jkl  1 k 


and 

/ C. ..  ,3  .u.a,uL  = 0 . 

Jv  ljkl  j i 1 k 

A prudent  look  at  the  definition  of  generalized  gradient  indicates 
the  symmetry  of  the  generalized  gradient;  therefor^  since  c is  positive 
definite,  3u  = 0 almost  everywhere,  and  by  n.j  = Cjjj-l3luk  we  have 
o * 0 a.e.  Since  3u  = 0 a.e.,  u is  unique  modulo  a rigid  displacement. 

Among  the  types  of  singularities  admissible  by  the  preceding 
theorem  are  any  stress  fields  in  t (ft)  (this  includes  the  crack  problem 


94 

which  has  stresses  of  order  0(r  ' ) as  r -*■  0),  however,  ueW^ii)  pre- 
cludes finite  discontinuities,  indeed  Smirnov  [75]  shows  if  ueW^si), 
then  u is  absolutely  continuous  almost  everywhere  in  ft.  Obviously, 
absolutely  continuous  implies  continuity. 

If  M is  a positive  measure  function  on  (X,M),  then  $:D  -*■  Y, 
where  (Y,||  • ||  ) is  a normed  linear  space  and  DcM,  is  a set  function 
on  D.  We  say  <p  is  absolutely  continuous  with  respect  to  p,  provided 
given  any  e > 0 there  exists  6 > 0 such  that 

U (E)|  < e whenever  u(E)  < 6,  ErD. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

This  work  has  produced  three  main  results— failure  zone  effects 
on  isochromatic  lines,  feasibility  of  laser  interferometric  techniques 
in  the  measurement  of  displacement  fields  near  crack  tips  in  solid 
propellant,  and  uniqueness  of  solution  to  a broad  class  of  elastostatic 
problems. 

Studies  of  the  isochromatic  fringe  patterns  in  an  infinite  sheet 
with  a semi-infinite  crack  under  equal  biaxial  loadings  indicate  that 
failure  zones  induce  tilting  of  isochromatic  fringe  patterns.  Addi- 
tional problems,  such  as  a finite  sheet  with  a small,  centered  crack 
(small  relative  to  sheet  dimensions)  under  equal  biaxial  loadings, 
need  to  be  solved  (numerically)  to  determine  the  extent  of  failure 
zone  effects. 

The  experimental  results  observed  and  recorded  herein  closely 
agree  (6%-10*)  with  theoretical  results  from  linear  theories  and  par- 
tially confirm  results  observed  by  others  [8].  There  were,  however, 
apparent  nonl inearities,  or  at  least  observations  not  explainable 
through  linear  theory.  To  further  substantiate  these  observations, 
additional  experimental  data  should  be  obtained  in  a more  controlled 
environment.  A loading  frame  that  would  eliminate  out-of-plane  and 
in-plane  translational  displacements  is  desirable.  Different  loadings 
on  the  sample  should  be  scrutinized,  e.g.,  the  combined  Mode  I and 
Mode  II  problem,  and  compared  to  analytical  solutions. 

The  uniqueness  of  solution  to  the  traction  boundary  value  problem 
proved  here  applies  to  a variety  of  traction  boundary  value  problems. 
Many  authors  [76-78]  treat  elliptic  boundary  value  problems  using  a 
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modern  approach.  (Uniqueness  for  the  generalized  Dirichlet  problem 
is  given,  for  example,  in  [76].)  It  is  the  way  the  boundary  conditions 
P (p.  91)  of  the  uniqueness  theorem  (as  opposed  to  normal  derivatives 
which  are  specified  in  the  general ;zed  Dirichlet  problem)  that 
separate  this  treatment  from  previous  work  [76-78]. 

In  summary,  it  is  the  view  of  this  author  that:  1)  further 
study  determining  the  extent  of  failure  zone  effects  should  be  con- 
ducted, 2)  further  refinements  are  needed  in  the  interferometric 
technique,  whose  feasibility  has  been  established  in  the  present 
work,  and  3)  uniqueness  of  the  displacement  boundary  value  problem 
should  be  investigated. 
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APPENDIX  A 

The  Equivalence  of  the  Singular  Terms  of  Williams  and  Irwin 


MODE  I 


From  equation  (23) 

2 2 

o =o  cos  o + o„  sin  o - 2cj  n sin  o cos  o 
x r e ro 


Equation  (7)  implies 

°r  ' 4^«['  5 C0$  ? * C°S  F °] 

if,  0 3 "1 

■ 3 cos  2 - cos  £ 0 I 


°0  = 


Tr0=  4r ' 


*i  r . e . 3 i 

T/2  " sin  2 " sin  2 0 
Equations  (A.l)  and  (A. 2)  give 


5 cos  cos2  o + cos  | 0 cos2  o 


J«[' 

0.2  3 2 

- 3 cos  2 sin  o - cos  ^ o sin  0 

0 3 

+ 2 sin  0 cos  o (sin  ^ + s1n  f °) 


Using 


cos  j 0 = cos  o cos  | - 2 sin2  cos  | 


0 O n 

sin^  + sin  | o = 2 sin  o cos  ^ 
sin2  | = 1(1  - cos  o)  , 


equation  (A. 3)  and  (A. 4)  give 


(A.l) 


(A. 2) 


(A. 3) 


( A.4) 


u 


1 

II 


103 


a1  o , o . 3 „ 

°x  = 172  cos  2 ” sin  2 sin  2 0 


(A. 5) 


Recalling  equation  (1), 


°x  ' (2,7)172 


o r i j o .3 

cos  2 1 - sin  2 sin  2 


(A. 6) 


Equations  (A. 5)  and  (A. 6)  imply 


51  = ' /In 


(A. 7) 


Equation  (A. 7)  could  be  found  by  using  the  expression  for  or 
t from  equation  (23). 


MODE  II 


From  Williams  (13] , 


- bl  f- 

' ^7?  L 

°r  ■ ^ [' 
t™  ■ 7/?  [■ 


5sin^  + 3sin|o 

2 c j 


bl  I"  o • 0 ^ . 3 „ 

o = — wo  " 3 sin  5 - 3 sin  w o 
r 4r’'Z  L ^ ^ 


(A. 8) 


3 sin  2'  “ 3 sin  2 A 


Equations  (A.l)  and  (A. 8)  give 
b,  r 


V72  s*n  2 cos^  0 + s’n  f 11  cos^  0 (A* 

■ 6 sin'*  0 cos  - 4 sin  ^ cos  ^ cos  0 (cos  + cos  ^ o)J 


(A. 9) 


where 


sin  I 0 = sin  I cos  0 + sin  0 cos  | . 


(A. 10) 


Equations  (A. 4)  (A. 9)  (A. 10)  imply 
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°x  = * rl/2  (si 


" f>  [ 


2 + COS  | COS  ^ 0 


]. 


Recalling  equation  (1), 
K 


II  o f,  , o 3 1 

0 _ sin  ^ 2 + cos  x cos  o o 

x (2vr)]/Z  2 l 2 2 J. 


Equations  (A. 11),  (A. 12)  imply 

b - Kn 
1 ' (2*)1'2  • 


(A. 11) 


(A. 12) 


(A. 13) 
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APPENDIX  B 

Linking  the  Failure  Stress  Distribution  with  TEXGAP 


The  following  material  properties  will  be  used  in  this  analysis, 

E (v ) . . .Young's  (shear)  modulus  , . 

I] 

v . . . .Poisson's  ratio. 

The  constants  in  (B.l)  are  related  through  [18] 


3 InT  • (B.2) 

Using  dimensional  analysis  for  a quarter  biaxial  strip  of  width  w. 


height  h,  crack  length  a,  and  grip  displacement  v^, 
v 


Kj  = Xu 


* 00 


(B.3) 


where  x is  a constant  of  proportionality. 

From  equations  (25a)  and  (25b)  corresponding  to  the  constant  and 
ramp  failure  stresses, 


= 2(2  a/n)1^2 

(constant  stress) 

(B.4) 

_ 4 , „ . ,1/2  n 

“3(2o/i,)  °m 

(ramp  stress) 

For  plane  stress  [79] 

G = K,2/(3u)  , (B. 5) 

where  G is  the  energy  release  rate  and  may  be  calculated  through 


(30] 


1 

2 


G > ~ AW/Aa  , (R.f) 

and  W is  the  strain  energy' per  unit  thickness,  and  is  obtained  from 

TEXGAP.  Equations  (B.3)  and  (B.5)  combine  to  give 

X2  vj2  h2  (B.7) 

G = ~ M 
3 a w 


I 


I 

] 

1 

I 

T 

1 

] 

I 
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2 v 2 2 

where  \ 9 h is  constant  for  a given  problem. 

3 a 

We  seek  an  expression  for  o . Equations  (R.3),  (B.4)  give 

A-l'-JUl  = 2 (2  a / it  )1/2  om 

w/ a 

v 1/2 

A n q h = 4 (2a\'  0 (B.8) 

w/a  3(ttf  m 

now  equations  (B.7),  (B.8)  combine  to  give 

om=  /{3Gvnv/(8n)}  (B.9) 

for  the  constant  failure  zone  stress  distribution  and 

om=  J {27Guit/(32  a ) > (B.10) 

for  the  ramp  failure  zone  stress  distribution. 
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APPENDIX  C 

The  Extended  Divergence  Theorem 

Let  n be  a closed  region  in  whose  boundary,  3 li , is  a rectifi- 
able Jordan  curve  and  let 

o - {(xj>  x^»  x3)  | y ( x ^ , Xg  )j5.  *3  _5.  $ (^i«  , ( x.| , Xp)  c n } , 

where  4>eC  1 ( n ) . Further  assume  4^  = 4>  i f and  only  if  (x^,  Xp)  c an 


Let  Sj  and  Sp  be  the  trace  of  the  surfaces  described  res- 
pectively through 


S^x)  = x1e1  + x2e2  + ^(x^  Xp)e3 
s2(x)  = x]e1  + x2e2  + ip ( x 1 9 x2)e3 
and  let  S = SjUSg. 


for  all  (x^ , x2 ) en 


Theorem.  (The  Extended  Divergence  Theorem).  Let  n be  a bounded 
regular  region  (hence  a Lebesque  measurable  set),  and  <(>,  u . o a 
scalar,  vector,  and  second-order  tensor  field,  differentiable  almost 


everywhere  on  ?}  with  derivatives  in  l (ii),  then  the  results  of  the 
classical  divergence  theorem  hold  with  the  proper  interpretations  of 
the  surface  integral . 


Proof 

We  need  only  consider  as  defined  in  the  preceeding  paragraph, 
and  will  only  show  the  proof  for  the  two-tensor  field  o as  different 
order  tensor  fields  are  proved  in  a similar  fashion. 

Now  « is  a Lebesgue  measurable  subset  of  R ; therefore,  by 
restricting v>  , the  Lebesgue  measure  on  R^,  to  the  measurable  sub- 
sets of  n,  a new  measure  space  is  obtained  in  the  obvious  fashion; 
indeed,  we  may  think  of  n as  a product  space  and  apply  fubini's 
Theorem  to  obtain 


*(xr  x2) 

J °i 3 3dl'  =J  1°. i3,3  dxdo 


% X?) 

where  X is  the  Lebesgue  measure  on  the  line  segment  between 


$(Xjt  Xp)  and  , x?)  and  o the  Lebesgue  measure  on  n.  By 

1 

oi3>3  exists  and  is  L (u)  therefore 

13,3  d"  = j [oi3(x1,  Xp,  ^(x-j,  x 2 ) ) - o^3(x^,Xp  '/'(x-j,  Xp))]do 

by  the  smoothness  of  Sj  we  obtain 

fa 13,3du  af^°  13  ° §1  ^ I JS,  I ‘ (°13  ° do' 

*'n  ' 2 

Define  the  vector  <1,3  (D,  0,oi3  ) then  by  the  convention 

established  in  the  uniqueness  chapter. 
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ON  THE  MECHANICAL  BEHAVIOR  OF  SOLID  PROPELLANT 
UNDER  TRANSIENT  TEMPERATURES 
R.  A.  Schapery 

ABSTRACT 

The  effect  of  simultaneous  cooling  and  straining  on  composite 
solid  propellant  Is  studied  experimentally  and  theoretically.  It  Is 
shown  that  the  use  of  isothermal  viscoelastic  relaxation  data  and  the 
linear  viscoelastic  theory  for  thermorheological ly  simple  materials 
results  in  stresses  which  greatly  underpredict  the  actual  stresses  in 
propellant  under  this  Important  type  of  loading.  An  improved  constitu- 
tive equation  for  predicting  the  mechanical  response  of  propellant  is 
then  proposed.  It  is  based  on  theoretical  results  obtained  from  non- 
equilibrium  thermodynamics  and  a viscoelastic  theory  of  microcrack 
growth.  Some  aspects  of  the  calculated  response  are  described  and 
shown  to  be  in  qualitative  agreement  with  experimental  data.  However, 
additional  experimental  results  are  needed  In  order  to  assess  the  ac- 
curacy and  range  of  validity  of  the  theory. 
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Composite  solid  propellant  consists  of  a viscoelastic  rubber  matrix 
which  Is  highly  filled  with  a broad  size  distribution  of  hard  particles 
of  ammonium  perchlorate  and  aluminum;  the  volume  percent  filler  is  ty- 
pically 75*.  The  weak  rubber  matrix  together  with  high  internal  stress 
concentrations  produced  by  the  filler  are  largely  responsible  for  caus- 
ing highly  nonlinear  behavior  over  practical ly  the  entire  range  of  ser- 
vice-induced stresses.  Fortunately,  there  now  exist  verified  engineer- 
ing methods  of  characterizing  and  stress  analyzing  propellant  in  terms 
of  "effective"  linear  properties  for  many  types  of  loading  [1].  Also, 
linear  viscoelastic  fracture  mechanics  has  met  with  considerable  suc- 
cess when  the  correct  stresses  are  employed  [2].  However,  it  is  not 
yet  possible  to  predict  the  correct  stresses  in  rocket  motors  subjected 
to  a large  temperature  drop  with  any  reasonable  degree  of  confidence. 

A few  ad  hoc  methods  exist  for  predicting  stresses  during  cooling,  but 
their  accuracy  is  highly  dependent  on  the  strain  and  temperature  his- 
tory and  the  particular  propellant  investigated  [1]. 

The  magnitude  of  the  discrepancy  between  linear  theory  and  actual 
nonlinear  nonisothermal  behavior  is  not  well -documented  in  the  litera- 
ture,and  therefore  we  discuss  this  aspect  of  the  problem  in  the  next 
section.  A generalized  constitutive  model,  which  is  believed  to 
correctly  account  for  the  main  physical  mechanisms  of  deformation, 
is  then  developed  and  discussed  in  the  last  section. 
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EXPERIMENTAL  RESULTS 

A series  of  biaxial  tensile  tests  on  1/8"  thick  strip  specimens, 
which  are  shown  in  Fig.  1,  and  uniaxial  tensile  tests  was  conducted 


t 


Figure  1.  The  strip-biaxial  sample. 

under  the  direction  of  Francis  [3].  The  specimens  were  subjected  to 
the  constant  strain  and  cooling  rates  shown  in  Fig.  2;  the  rates  were 
applied  simultaneously  in  order  to  simulate  the  type  of  loading  which 
is  experienced  by  propellant  when  in  a case-bonded  motor  under  cooling 
conditions.  The  abscissa  in  Fig.  2 is  a non-dimensional  temperature 
defined  as 


where  Tq  is  the  initial  temperature  (=  70°F),  T is  the  instantaneous 
temperature  and  T is  a constant  (“  -145°F);  AT  is  a parameter  which 
is  employed  iri  the  linear  theory.  The  ordinate  is  the  ratio  of  the  ex- 
perimentally derived  engineering  stress,  a,  to  the  stress,  Oy,  pre- 
dicted from  linear  theory  for  a thermorheologically  simple  material. 
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Ratio  of  experimental  to  predicted  stresses  in  HTPB  solid 
propellant  under  simultaneous  cooling  and  straining. 
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Also  shown  In  Fig.  2 Is  this  stress  ratio  at  the  case  bond  line;  the 
stress  o Is  an  average  quantity  obtained  from  stress  gages  In  small 
steel-cased  motors  (4.5  Inch  outside  diameter  x 24  Inches  long  with  ap- 
proximately 1 Inch  diameter  circular  port). 

The  stress  Oj  was  calculated  using  the  method  In  Section  5. 5. 3. 3 of 
[1].  This  method  employs  the  uniaxial  relaxation  modulus  E ■ £(4)  and 
time- temperature  shift  factor  a^  = a-j-(T);  4 = t/a^  Is  "reduced  time". 

The  modulus  and  shift  factor  are  shown  In  Fig.  3 and  Fig.  4 respective- 
ly; the  crack  growth  data  cited  In  Fig.  4 are  discussed  In  [2]. 

It  Is  Indicated  In  Fig.  3 that  the  relaxation  modulus  was  obtained 
at  a strain  of  3%.  Inasmuch  as  variable  straining  tests  [3]  revealed 
that  the  material  is  nonlinear  at  this  strain  level  (by  checking  against 
the  superposition  principle),  two  other  types  of  tests  were  conducted 
at  Texas  A&M  to  obtain  the  relaxation  modulus.  The  duomorph  test  [4] 
consisted  of  placing  a thin  piezoelectric  sandwich  disc  against  a block 
of  propellant  and  noting  the  effect  of  the  propellant  on  vibrational 
bending  strain;  linear  viscoelastic  analysis  yielded  the  relaxation 
modulus  [4],  The  other  type  of  test  consisted  of  resting  long  circular 
cylinders  under  constant  load  on  the  same  block  of  material  and  observ- 
ing the  amount  of  time- dependent  Indentation;  again,  linear  analysis 
yielded  the  desired  modulus  [5].  The  average  strains  In  the  propellant 
for  the  duomorph  and  Indentation  tests  were  approximately  0.0H  and  3*, 
respectively.  The  data  from  all  three  tests  are  In  excellent  agreement, 
as  shown  In  Fig.  3.  That  the  duomorph  data  were  obtained  at  very  small 
values  of  time  at  room  temperature,  and  the  fact  that  they  agree  out  to 
the  low  temperature  end  (»  -65°F)  of  the  relaxation  modulus,  serves  as 
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a direct  check  on  the  validity  of  the  shift  factor  in  Fig.  4.  This 
agreement  between  the  results  of  all  three  tests  is  very  surprising; 
but  It  can  be  "explained"  by  using  the  constitutive  equation  discussed 
in  the  next  section. 

Returning  to  Fig.  2,  it  is  seen  that  the  experimentally  obtained 
stress  greatly  exceeds  the  predictions.  There  appears  to  be  consider- 
able specimen  - to  - specimen  property  variation  (cf.  the  first  two 
tests  In  the  table),  but  it  Is  clearly  not  sufficient  to  account  for  the 
discrepancy  between  theory  (<jy)  and  experiment  (a). 

We  shall  discuss  next  a constitutive  equation  which  Is  believed  to 
be  sufficient  for  removing  this  difference. 
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A CONSTITUTIVE  EQUATION  FOR  TRANSIENT  TEMPERATURES 


Let  us  first  consider  a linear,  uniaxial  stress-strain  equation  for 
a special  type  of  thermo rheological ly  complex  behavior  [6],  which  was 
originally  derived  from  thermodynamic  principles: 
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ar./EU  - S') 


de 
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(2) 


where 


5 = £(t) 


i = C(t) 


(3) 


(Note  that  the  reduced  time  s reduces  to  that  used  In  Fig.  3 for 
constant  temperature.)  Also  eo  is  the  total  strain  less  that  due  to 
thermal  expansion. 


ea  = c - a(T  - T0)  (4) 

where  a is  the  expansion  coefficient,  and  aQ  * aQ(T)  Is  a temperature- 
dependent  material  function.  If  one  were  to  use  a mechanical  model 
consisting  of  springs  and  dashpots  to  represent  Eq.  (2),  It  can  be 
easily  shown  that  all  of  the  spring  moduli  would  be  proportional  to  aG 
and  all  of  the  dashpot  viscosities  would  be  proportional  to  the  product 
aTaG.  If  aQ  * 1 these  moduli  are  independent  of  temperature,  and  Eq. 
(2)  reduces  to  the  constitutive  equation  for  a thermorheologlcally  sim- 
ple material. 

For  an  isothermal  relaxation  test,  Eq.  (2)  becomes 
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o * a E({)e  (5) 

b ° 

where  £ = t/a-r.  If  we  assume  that 

E(0  - E^""  (6) 

where  E1  and  n are  constants,  then  Eq.  (5)  can  be  written  as 

0 • E (£)e0  (7) 

where 

1 = t/(aTaG1/n)  (8) 

Thus,  as  observed  previously  [6],  if  one  obtains  experimental  values  of 
the  relaxation  modulus,  o/eo,  at  different  temperatures,  the  standard 
procedure  of  shifting  data  to  obtain  a master  curve  will  yield  only  a 
composite  shift  factor,  ac=  aTaG1/n;  thus,  for  power  law  materials, 
isothermal  data  are  not  sufficient  to  evaluate  all  material  functions  in 
Eq.  (2).  On  the  other  hand  for  many  transient  temperature  conditions, 
such  as  a constant  cooling  rate,  one  finds  that  the  stress,  Eq.  (2),  is 
sensitive  to  the  values  of  aG  and  aT  separately;  therefore,  nonisother- 
mal  testing  may  be  employed  to  complete  the  characterization. 

These  observations  are  highly  relevant  to  solid  propellant  since 
its  modulus  can  be  approximated  by  a power  law  over  a wide  range  of 
reduced  time;  e.g.,  the  slope  of  the  log-log  plot  In  Fig.  3 is  a very 
weak  function  of  5.  Furthermore,  a study  of  Eq.  (2)  reveals  that  It  is 
possible  to  select  functions  aG  for  which  o/oj  » 1 by  letting  ag 
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Increase  with  decreasing  temperatures;  this  can  be  achieved  without 
difficulty  without  affecting  the  Isothermal  behavior.  We  conclude, 
therefore,  that  at  least  a portion  of  the  discrepancy  between  experi- 
ment and  the  conmonly  used  theory  (l.e.  aQ  ■ 1 ) may  be  due  to  a type  of 
thermorheologlcal  complexity  for  which  elastic  moduli  increase  as  the 
material  Is  cooled. 

This  type  of  generalization  Is  apparently  not  sufficient  since  the 
stress  ratio  In  Fig.  2 varies  with  the  straining  and  cooling  rates;  some 
of  the  differences  between  the  curves  Is  due  to  material  variability,  as 
discussed  above,  but  probably  not  all.  Also,  as  noted  previously.  Iso- 
thermal tests  reveal  that  the  material  is  nonlinear. 

A significant  If  not  major  portion  of  this  nonlinear  behavior  may 
be  due  to  microcracking  [7,8],  Here  we  shall  draw  upon  the  model  devel- 
oped in  [9]  In  order  to  Incorporate  this  phenomenon  in  the  nonisothermal 
constitutive  theory.  From  [9],  as  generalized  by  means  of  Eq.  (2), 
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where  D(c)  Is  the  uniaxial  creep  compliance, and  for  a non-aging  material 
F Is  a function  of  the  "damage"  parameter 
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where  q is  a positive  constant.  The  microstructure  serves  to  define  F 
through  the  distribution  of  Initial  stress  Intensity  factors  for  the 
microflaws  and  their  softening  effect  during  growth;  however.  In  prac- 
tice, F can  be  found  from  experimental  stress-strain  data.  It  Is  shown 
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In  [9]  that  dF/dSQ  > 0 in  view  of  Its  relation  to  the  microcracking. 
Each  microcrack  was  assumed  to  obey  the  growth  rate  equation. 


t ■ V on 

where  A and  q are  positive  constants,  and  Kj  is  the  opening  mode  stress 
intensity  factor. 

Equation  (9)  may  be  inverted  exactly  to  find 
ar  r de 

°--t/eU  - *')  --37  cIt  (12) 

Thus,  we  see  that  Eq.  (2)  is  contained  in  nonlinear  Eq.  (12). 

One  further  generalization  of  Eq.  (2)  will  be  made  to  account  for  the 
effect  of  strain  on  Eq.  (11).  Specifically,  Eq.  (11)  was  found  to  ap- 
ply to  macrocracks  if  a retarding  effect  of  tensile  strain  is  incorporat- 
ed in  the  form  [2], 


£-=A(Kj/f)q  (13) 

where  f * f(eo)  is  an  increasing  function  of  overall  specimen  strain; 
for  a large  amount  of  data  it  was  found  that 

C e ® , 0 < < 0.085 

f » 0 0 (14) 

constant,  eo  > 0.085 

where  s « 0.7.  It  will  be  assumed  that  strain  affects  microcrack  growth 
in  the  same  way,  and  therefore  f must  be  included  in  the  constitutive 
theory.  A review  of  the  derivation  in  [9]  shows  that  the  only  modifica- 
tion needed  is  the  replacement  of  SQ  by  the  more  general  damage  para- 
meter. 
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So  far  we  have  discussed  Eq.  (9)  as  If  It  were  limited  to  uniaxial 
tensile  loading  of  a specimen  containing  only  opening-mode  cracks.  How- 
ever, our  results  are  actually  more  general  than  this.  Specifically,  if 
proportional  loading  exists  and  if  mixed-mode  crack  growth  obeys  an 
equation  of  the  same  form  as  Eq.  (13),  but  with  Kj^  replaced  by  any 
homogeneous  function  of  the  three  stress- intensity  factors,  then  Eq.  (2) 
still  applies.  In  this  more  general  case,  eQ  and  o may  be  intepreted  as 
the  maximum  principal  strain  and  stress,  respectively,  and  EU)  is  the 
linear  viscoelastic  relaxation  modulus  for  the  particular  ratio  of  prin- 
cipal stresses  under  consideration.  (This  interpretation  will  be  used  in 
the  rest  of  the  analysis.)  The  function  F is  also  dependent  on  the 
stress  ratios;  but,  for  given  ratios,  F = F (S)  and  dF/dS  > 0.  Pro- 
portional loading  may  not  be  exactly  satisfied  for  the  strip  geometry  in 
Fig.  1,  but  because  the  lateral  strain  vanishes,  the  results  apply  to 
this  problem.  Of  course,  there  is  no  basis  for  assuming  F is  the  same 
for  uniaxially  stressed  bar  and  strip  specimens. 

Let  us  now  rewrite  Eq.(12)  in  a more  useful  form  by  solving  for 
stress  in  terms  of  strain  and  temperature  history.  First,  define 


(16) 


which  is  the  stress  for  a linear  viscoelastic,  thermorheologically  sim- 
ple material.  Also,  define  the  ratio. 
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Equation  (23)  Is  the  desired  result  as  It  expresses  the  stress  In 
terms  of  the  strain  and  temperature  histories;  these  histories  appear 
In  Eqs.  (16)  and  (20).  The  function  F^S^)  Is  directly  related  to  the 
initial  microstructural  state  of  the  material;  It  is  to  be  found  by 
applying  Eq.  (23)  to  experimental  data  rather  than  carrying  out  detailed 
calculations  from  microstructural  Information.  We  should  add  that  on 
the  basis  of  some  existing  data  the  function  f was  assumed  to  be  a func- 
tion of  applied  strain  but  not  of  temperature  [2],  If,  however,  the 
shift  factor  for  microcracks  is  different  from  that  for  overall  mechan- 
ical behavior  and/or  If  a function  analogous  to  as  Is  used  in  the  pre- 
diction of  crack  speed  from  creep  compliance  Information  [9],  then  f 
will  depend  on  temperature.  Consequently,  for  the  sake  of  generality, 
we  shall  suppose  f * f(eo,T). 

Let  us  now  discuss  the  previously  reported  experimental  data  in 
light  of  Eq.  (23).  First  Introduce  the  secant  moduli. 
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and  obtain 
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Now,  If  f ~e  , Eqs.  (20)  and  (25)  show  that  the  nonlinear  secant 
0 

modulus  E Is  a function  of  the  history  of  E.  but  Is  Independent  of  the 
n 1 

strain  level.  Thus,  for  example,  the  relaxation  modulus  Is  predicted 
to  be  Independent  of  strain  even  though  the  material  Is  nonlinear;  this 
behavior  often  exists  for  solid  propellant  at  small  strains.  Also,  for 
a constant  amount  softening  due  to  damage  (l.e.  * constant). 
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Eq.  (25)  shows  that  linear  viscoelastic  modulus  Interrelationships  apply 
to  the  nonlinear  modulus,  which  was  found  to  be  true  for  the  HTPB  propel- 
lant (cf.  Fig.  3).  This  type  of  behavior  will  exist,  at  least  as  an 
approximation  ,1f  there  is  a range  for  which  Fi  is  a weak  function  of 
damage  (i.e.,  a range  for  which  there  Is  a negligible  or  only  a small 
change  with  the  number  and  size  of  microcracks).  For  this  situation  the 
ratio  of  the  measured  relaxation  modulus  to  the  linear  viscoelastic 
modulus  is  constant.  The  modulus  in  Fig.  3, for  example,  would  then 
be  a constant  fraction  of  the  true  linear  modulus.  Inasmuch  as 
ay(cf . Fig.  2)  was  calculated  using  the  data  in  Fig.  3,  we  conclude  that 
o^/oj  = K,  where  K > 1 is  a constant.  The  theory,  Eq.  (23),  correspond- 
ing to  Fig.  2 then  is  written 


If  the  softening  due  to  damage  were  the  same  for  all  data  in  both  Figs. 

2 and  3,  we  would  have  o/oj  = aQ;  the  spread  in  the  data  in  Fig.  2 would 
then  be  due  solely  to  material  variability,  and  all  dependence  on  temper- 
ature could  be  attributed  to  a^.  On  the  other  hand,  F^  is  probably  not 
constant  and,  in  fact.  Is  believed  to  be  the  cause  of  the  decrease  in 
some  of  the  data  in  Fig.  2 at  large  values  of  aT^.  Most  of  the  propel- 
lant In  the  subscale  motor  Is  at  relatively  very  low  distortion,  and 
for  this  case  S should  exhibit  much  less  of  an  increase  than  for  the 
strip  tests;  this  observation  Is  consistent  with  the  data  in  Fig.  2. 

Another  interesting  phenomenon  is  observed  when  a tensile  bar  of 
solid  propellant  Is  held  at  a constant,  positive  strain  for  a period 
of  time  and  then  cooled  at  a constant  rate:  the  rate  of  stress 
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increase  is  often  much  greater  than  when  the  strain  is  zero  [3].  This 
behavior  is  predicted  by  Eq.  (23)  if  the  factor  aG  increases  significant- 
ly with  decreasing  temperature.  Furthermore,  the  data  from  such  tests 
are  very  useful  for  obtaining  this  material  function.  Indeed,  various 
data  are  now  being  generated  by  Francis  [3]  in  order  to  establish  both 
functions  aG  and  F^  and  check  the  predictions  of  the  constitutive  equa- 
tion (23). 
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SOME  VISCOELASTIC  CRACK  GROWTH  RELATIONS  FOR  ORTHOTROPIC 


AND  PRESTRAINED  MEDIA* 
by  G.  S.  Brockway  and  R.  A.  Schapery 

ABSTRACT 

Equations  are  developed  for  predicting  opening  and  sliding  modes  of 
crack  growth  along  planes  of  geometric  symmetry  in  viscoelastic  ortho- 
tropic media  with  and  without  large  prestrains.  Except  for  the  small 
failure  zone  at  the  crack-tip,  the  material  is  assumed  to  be  linearly 
viscoelastic  with  respect  to  the  changes  in  stress  and  strain  which  occur 
during  crack  propagation.  For  an  orthotropic  body  in  plane  strain,  an 
observation  of  Biot  relating  the  Fourier  transformed  tractions  and  dis- 
placements on  the  surface  of  any  homogeneous,  elastic  half-space  is  used 
in  conjunction  with  Graham's  extended  correspondence  principle  to  gene- 
rate the  viscoelastic  crack  face  displacements. Appl ication  of  a suitable 
energy  criterion  for  failure  leads  to  a nonlinear  integro-differential 
equation  for  the  crack  tip  speed  in  terms  of  material  and  loading  para- 
meters. The  resulting  equation  is  of  the  same  generic  form  as  in  the 
previously  published  isotropic  case  except  that  the  isotropic  creep  com- 
pliance is  replaced  by  effective  compliances  formed  from  orthotropic 
moduli.  Thus,  just  as  in  the  Isotropic  case,  a much  simpler,  approximate 
crack  growth  rate  equation  is  deduced  from  the  assumption  these  effective 

♦Published  originally  as  Texas  A&M  Univ.  Rept.  No.  MM- 3064-75-3  (AFOSR 
TR-76-1148)  1975;  to  be  published  In  Engineering  Fracture  Mechanics, 1978. 


127 


compliances  have  small  curvature  when  plotted  logarithmically  against 
logarithmic  time.  Included  is  a numerical  sample  using  constitutive 
properties  of  a fiber- reinforced  plastic.  Extension  to  fracture  of 
prestrained  media  and  crack  growth  between  certain  types  of  dissimilar 
media  is  then  made.  Again  the  growth  equations  have  the  same  basic 
form,  but  for  a prestrained  medium  the  effective  compliances  involve  the 
incremental  moduli.  The  results  indicate  that  crack  growth  is  acceler- 
ated for  initial  compressive  stress  states  that  create  certain  surface 
instabilities  ahead  of  the  crack  tip.  Finally,  it  is  argued  that  these 
crack  velocity  relations  can  be  employed  in  plane  stress  problems  and  in 
many  three-dimensional  cases,  and,  with  a small  change,  can  be  used  to 
predict  the  time  of  fracture  initiation. 

CONCLUSIONS 

Equations  for  predicting  crack  growth  in  opening  and  sliding  modes 
in  orthotropic  media  have  been  developed.  Although  the  emphasis  was  on 
homogeneous,  orthotropic  media  in  plane  strain,  It  was  shown  that  for 
the  opening  mode  these  crack  growth  relations  carry  over  without  change 
to  three-dimensional  problems  If  the  material  is  Isotropic  in  the  crack 
plane.  The  relations  were  also  extended  to  problems  involving  crack 
propagation  between  certain  types  of  dissimilar  orthotropic  media. 

Furthermore,  by  means  of  an  argument  analogous  to  that  used  pre- 
viously for  Isotropic  media,  one  can  show  that  the  plane-strain  crack 
growth  relations  apply  even  when  the  global  state  of  the  body  Is  one  of 
plane  stress;  it  is  only  necessary  that  the  failure  zone  be  small 
enough  for  the  neighborhood  of  the  crack-tip  to  be  In  plane  strain.  The 


assumption  of  crack  plane  Isotropy  Is  not  needed. 

Except  for  a small  modification.  It  was  demonstrated  that  the  crack 
growth  relations  in  terms  of  Incremental  properties  apply  to  homogeneous- 
ly prestrained  media  If  linearity  exists  with  respect  to  the  incremental 
stresses  and  strains.  Although  no  theoretical  limits  were  placed  on  the 
magnitudes  of  the  initial  strains,  we  did  assume  that  the  Initial  maximum 
shear  stress  was,  at  most,  of  the  same  order  of  magnitude  as  the  incre- 
mental stresses  associated  with  crack  propagation.  Allowing  for  large 
initial  strains,  while  limiting  the  initial  shear  stress  and  Incremental 
strains,  are  not  necessarily  incompatible  conditions  with  viscoelastic 
materials  which  undergo  cracking.  Consider,  for  example,  the  situation 
in  which  the  body  is  strained  very  slowly  and  assume  that  crack  growth 
does  not  occur  until  the  strains  are  quite  large  and  possibly  considerable 
material  nonlinearity  is  exhibited.  Then,  suppose  that  the  crack  grows 
during  an  interval  of  time  which  Is  small  relative  to  the  time  of  applica- 
tion of  the  externally  applied  strains.  The  incremental  strains  associ- 
ated with  the  unloading  process  of  crack  propagation  will  be  relatively 
small  (except  possibly  for  the  material  very  close  to  the  crack-tip)  If 
these  two  time  scales  are  far  apart  and  the  body  Is  highly  viscoelastic; 
as  just  one  application,  we  believe  these  requirements  are  quite  realistic 
in  certain  situations  encountered  with  solid  propellant  rocket  fuel. 
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A THEORY  OF  TIME-DEPENDENT  BONDING 
FOR  VISCOELASTIC  MEDIA* 

by  R.  A.  Schapery 

ABSTRACT 

A theory  is  developed  for  the  prediction  of  time-depending  bonding 
of  surfaces  of  similar  or  dissimilar  linear  viscoelastic  materials. 
Starting  with  an  extended  correspondence  principle,  viscoelastic  stresses 
and  displacements  in  the  neighborhood  of  a moving  contact  edge  are 
derived.  These  results  are  then  employed  in  a local  energy  balance 
which  leads  to  an  equation  for  predicting  the  rate  of  bonding  as  a func- 
tion of  the  local  displacement  intensity  factor.  A comparison  of  true 
and  apparent  energies  for  bonding  and  fracture  is  made.  Finally,  criteria 
are  given  for  determining  whether  or  not  bonding  can  occur,  and  then  the 
prediction  of  the  time  required  for  complete  bonding  or  healing  of  cracks 
is  examined. 

CONCLUSIONS 

A theory  of  time-depending  bonding  has  been  developed  which,  in 
certain  respects,  is  described  by  relations  which  are  similar  to  those 
of  time-dependent  fracture.  As  one  particular  result,  it  was  shown  that 
the  apparent  or  pseudo  fracture  energy,  which  is  dependent  on  the  speed 
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of  the  crack,  Is  typically  much  larger  than  the  speed-dependent  pseudo 
energy  of  bonding.  Additionally,  criteria  for  bonding  or  crack  healing 
were  derived  and  compared  with  the  conditions  under  which  crack  growth 
occurs. 

A simple  extended  correspondence  principle  was  employed  In  the 
development  of  the  theory  which  permitted  the  direct  use  of  elasticity 
solutions  In  predicting  viscoelastic  stresses  and  displacements.  Although 
this  principle  limits  us  to  certain  types  of  boundary  conditions.  It  is 
believed  the  theory  Is  applicable  to  a large  number  of  practical  problems 
dealing  with  mlcrccrack  healing  in  composites  and  adhesion  phenomena. 


I 
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A METHOD  FOR  PREDICTING  CRACK  GROWTH  IN 
NONHOMOGENEOUS  VISCOELASTIC  MEDIA* 

by  R.  A.  Schapery 

ABSTRACT 


Equations  are  developed  for  predicting  crack  growth  in  the  opening 
mode  for  quite  general  situations,  including  many  quasi-static  and  dyna- 
mic problems  involving  moisture  and  temperature  gradients  in  monolithic 
and  composite  materials.  Except  for  the  small  zone  of  failing  material 
at  the  crack  tip,  the  body  is  assumed  to  be  linearly  viscoelastic.  A 
method  for  obtaining  viscoelastic  stresses  and  displacements  from  elastic 
solutions  is  first  described.  The  traction  boundary  condition  for  the 
crack  faces  is  not  in  general  satisfied  by  these  results.  However,  it  is 
shown  that  by  modifying  the  failure  zone  this  condition  can  be  met.  An 
integral  equation  for  the  stress  in  the  modified  failure  zone  is  then 
derived.  Approximate  analysis  similar  to  that  used  previously  by  the 
author  in  stress  analysis  is  then  employed  to  solve  the  integral  equation 
and  develop  relatively  simple  equations  for  predicting  crack  speed;  these 
equations  relate  crack  speed  to  the  stress  intensity  factor  in  a suitably 
defined  elastic  body. 


♦Published  originally  In  Proceedings  of  Conference  - Environmental  Degra- 
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CONCLUSIONS 

The  problem  of  predicting  crack  speed  in  nonhomogeneous  viscoelastic 
bodies  has  been  reduced  to  quasi-elastic  analysis  through  a series  of  ap- 
proximations. The  primary  result  is  analogous  to  that  derived  previously 
for  situations  in  which  the  stress  intensity  factor  is  independent  of  com- 
pliance. For  composite  materials  and/or  when  the  viscoelastic  properties 
vary  throughout  the  body  (due  to,  for  example,  moisture  and  temperature) 
the  speed  controlling  compliance  Cy(ta)  is  that  for  the  material  in  the 
neighborhood  of  the  crack  tip.  Through  a proper  selection  of  this  com- 
pliance the  present  results  will  apply  to  cohesive  cracks  in  locally 
isotropic  media  and  to  cohesive  and  adhesive  cracks  in  certain  locally 
isotropic  and  orthotropic  media. 

The  far-field  viscoelastic  properties  affect  crack  growth  through 
the  quasi-elastic  stress  intensity  factors  Kje( t^ ) and  Kje(3ta);  these 
factors  are  the  instantaneous  values  for  an  elastic  material  having  the 
applied  loads  defined  in  the  paper  and  compliances  equal  to  viscoelastic 
creep  compliances  which  are  evaluated  at  the  "local"  times  t and  3tQ. 

The  functions  Cv  and  KjC  may  depend  on  actual  time  t in  addition 
to  local  time  ta,  although  the  former  time-dependence  is  not  explicitly 
shown  in  these  results.  This  dependence  on  actual  time  may,  for  example, 
be  due  to  changes  in  applied  loads,  temperature,  and  moisture  and/or 
changes  in  crack  size  and  speed.  However,  any  changes  which  do  occur 
should  be  small  during  the  time  required  for  the  crack  to  propagate  a 
distance  of  approximately  five  to  ten  times  the  length  of  the  failure 
zone,  a;  further  study  is  needed  to  determine  the  effect  of  relaxing 


this  condition. 
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Inertia  effects  associated  with  very  high  crack  speed  have  been 
neglected  in  the  development  of  explicit  results.  It  is  anticipated, 
however,  that  effects  of  high  speed  can  be  taken  into  account  using  the 
procedure  described  herein  because  the  dynamic  crack-tip  stress  field 
for  an  elastic  material  is  completely  characterized  by  a time-dependent 
stress  intensity  factor. 

The  present  analysis  enables  one  to  use  experimental  crack  growth 
data,  a = function  (Kje),  obtained  with  special  loading  histories  and 
uniform  temperature  and  moisture  conditions,  in  the  prediction  of  growth 
in  complex  structures  under  quite  general  loading  histories  and  nonuniform 
conditions.  However,  in  contrast  to  problems  for  which  Kje  does  not 
depend  on  compliances,  the  size  of  the  failure  zone  a is  explicitly 
needed  because  it  affects  the  value  of  the  stress  intensity  factor. 


